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Abstract 

It is shown how pair production itself would prevent the astrophysical for- 
mation of dyadospheres, hypothetical regions where the electric field exceeds 
^ ' the critical value for rapid Schwinger pair production. Pair production is a 

self-regulating process that would discharge a growing electric field, in the ex- 
ample of a hypothetical collapsing charged stellar core, before it reached 6% 
of the minimum dyadosphere value, keeping the pair production rate more 
than 26 orders of magnitude below the dyadosphere value, and keeping the 
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1 Introduction 



Ruffini and his collaborators 000111313111013110111111121113131131131111113 
C3I23I2I1E21E3I21E3E3EZ1EH] have proposed a model for explaining gamma 
ray bursts that presumes the initial existence of what they call a dy ado sphere, a 
macroscopic region of spacetime where the electric field exceeds the critical electric 
field E c = m 2 e c , /(Tie) ~ 1.32 x 10 16 V/cm for sufficiently rapid Schwinger pair 
production [29J. The difficulty of producing these large electric fields is a problem 
with this model that has not been adequately addressed. 

There are at least two strong reasons for doubting that such large electric fields 
can develop over astrophysical scales (i.e., over length scales much larger than the 
collision regions of individual charged particles). First, it would be very difficult 
to develop sufficient charge imbalance for macroscopic electric fields to produce 
significant numbers of pairs. Second, even if macroscopic pair production could 
somehow be achieved, I shall show in this paper that this process is sufficiently 
self-regulating that it prevents the electric field from achieving a value that would 
produce pairs at even 10~ 26 that of dyadosphere models. I conclude that it is highly 
implausible that dyadospheres can form in outer space, and therefore invoking them 
for models of gamma ray bursts is not viable. 

The first reason for being extremely doubtful of the existence of astrophysical 
dyadospheres is that it is very difficult for a large charge imbalance to develop 
astrophysically, because of the very high charge-to-mass ratio of elementary particles. 
For example, the ratio of the electrostatic repulsion to the gravitational attraction 
of two protons, each of charge q and mass m p , is the square of their charge-to-mass 
ratio in Planck units, which is 



This implies that if one had a spherical object, such as a stellar core, with a positive 
charge-to- mass ratio Q/M greater than the inverse of the charge-to- mass ratio of 
the proton, m p /q w 9 x 10~ 19 , the electrostatic repulsive force on the protons at 
the surface would be greater than the gravitational attractive force, so such protons 
would most likely be ejected. (If the object had a negative charge, electrons of mass 
m = m e would be expelled if —Q/M > m e /q « 4.9 x 10~ 22 , lower by the factor 
of the mass ratio of the proton and the electron, m p /m e m 1836, so the maximum 
value of the charge of such an object would be even less if it were negative.) 




(1) 
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If one takes the mass-to-charge ratio of the proton to be a rough estimate of the 
maximum charge-to-mass ratio of an astrophysical object (or else gravity would not 
be strong enough to hold in the protons that make up the excess charge) , then using 
the formulas of the succeeding sections, one can readily calculate that at the surface 
of a spherical object of radius R and mass M, the ratio of the electric field value, E, to 
the critical field value of a dyadosphere, E c = m 2 c 3 /(Tiq) rs 1.32328539 x 10 16 V/cm, 
is 

E_ J^m, / Mq\ /2GM\ 2 / M&\ 

E c ~ 4GM Q m e I M J U 2 /J { M ) ' {2) 

Therefore, if protons can be ejected from an astrophysical object whenever the 
electrostatic repulsion exceeds the gravitational attraction, then the electric field is 
constrained to be more than 13 orders of magnitude smaller than the critical value 
for a dyadosphere (if the mass is greater than 1.2 solar masses, which would be a 
conservative lower limit on any mass that could contract to 2GM/(c 2 R) ~ 1). For 
a negatively charged object, the corresponding limit would be more than 16 orders 
of magnitude smaller than the critical value for a dyadosphere. 

Although it seems very unlikely to work, one might seek to evade the electrostatic 
expulsion of protons by postulating that they are bound by nuclear forces to an 
astrophysical object, such as a collapsing neutron star core. The critical electric 
field E c that Ruffini and his collaborators [HI2l3l210000l3IIiBimiI2II31 
HH H3 IIHl HZl HH [HI |201 12H 1221 12S1 12H I23 123 12ZI I2HI used to define the minimum 
value for a dyadosphere would give the electrostatic force on an electron or proton 
of magnitude F c = qE c = (m e c 2 ) 2 / (he) ~ 0.00132 MeV/fm (2n times the rest mass 
energy m e c 2 of an electron divided by the Compton wavelength 2irhc/ (m e c 2 ) of 
the electron), whereas nuclear energies of the order of an MeV (e.g., the deuteron 
binding energy of about 2.225 MeV) over length scales of the order of a fermi (e.g., 
the proton charge radius of about 0.87 fm) would give nuclear forces of the order of 
a few MeV/fm, about 3 orders of magnitude larger than the electrostatic force of a 
minimal dyadosphere on a proton. Then one might suppose that the gravitational 
force can be provided by the neutrons of a collapsing neutron star core, and then a 
tiny fraction of the nucleons are in the form of protons that are bound by nuclear 
forces to the surface of the neutron star core, a sufficient fraction to give the electric 
field of a dyadosphere. 

My suspicion is that such a configuration would be highly unstable to pieces of 
the charged surface breaking off and being ejected by the huge electrostatic forces on 
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them, particularly as the core collapses and cannot maintain a fixed configuration 
of its nucleons (say a crystalline structure with strong nuclear forces between all 
neighboring nucleons). I would think that it is more reasonable that any strong 
nuclear crystalline structure would develop cracks that would make the effective 
nuclear forces too weak to prevent the charged surface of the core from breaking off 
and being ejected by an electric field anywhere near dyadosphere values. However, 
I do not have a firm result of when this would definitely happen, and therefore I 
was led to do the calculations below with the second mechanism (self-regulation) for 
preventing the occurrence of a dyadosphere. 

As an aside, let me describe one potential mechanism for ejecting the charge if 
protons are bound by nuclear forces to a collapsing core, though it proved not to be 
sufficient. This was the idea that radiation would convert the protons to neutrons, 
positrons, and anti-neutrinos, and then the positrons (not held to the core by nuclear 
forces) would immediately be ejected by the strong electrostatic forces. Presumably 
there must be enough radiation present to prevent the electrons in the inevitable 
plasma around the core from being pulled back to it to neutralize it. Then the hope 
was that radiation strong enough to keep away the electrons would also be strong 
enough to convert the charge carriers from protons (bound to the core) to positrons 
(not bound). However, this would be a weak interaction rate, and it turned out 
that the cross section [30J is apparently not large enough for the minimal amount 
of thermal radiation needed to keep the electrons at bay to be sufficient to convert 
the protons to positrons within the collapse time of the core. 

Therefore, I was led to do an analysis of the self-regulation of the pair production 
process itself, which, as we shall see below, will discharge any growing electric field 
well before it reaches dyadosphere values. This occurs essentially because astrophys- 
ical length scales are much greater than the electron Compton wavelength, which is 
the scale at which the pair production becomes significant at the critical electric field 
value for a dyadosphere. Therefore, the electric field will discharge astrophysically 
even when the pair production rate is very low on the scale of the electron Compton 
wavelength. 

The calculations below lead to the conclusion that it would be very difficult 
astrophysically to achieve, over a macroscopic region comparable to the size of a 
black hole or larger star, electric field values greater than a few percent of the 
minimum value for a dyadosphere, if that. The Schwinger pair production itself 
would then never exceed 10~ 26 times the minimum dyadosphere value. 
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2 Schwinger discharge of a charged collapsing core 



In this section we shall analyze the pair production and discharge of an electric field 
produced by the collapse of a hypothetical charged sphere or stellar core, ignoring 
the processes already discussed that would most probably cause almost all of the 
excess charge of the sphere to be electrostatically ejected. 

For a sphere collapsing in finite time, there is only a finite time for the electric 
field to be discharged by the Schwinger pair production process, so the discharge is 
never complete, but instead leaves a residual electric field at each moment of the 
collapse. We shall calculate an upper limit for this value and show that it is always 
more than a factor of 18 less than the minimum value for a dyadosphere. Because 
the pair-production rate is exponentially damped by the inverse of the electric field 
value, the pair production rate itself never exceeds a value that is more than 26 orders 
of magnitude below that of a dyadosphere. (The factor of the order of 10 26 , which 
we shall derive below, comes mainly from the fine structure constant multiplied by 
the square of the ratio of the Schwarzschild radius of a solar mass black hole to the 
Compton wavelength of an electron, which is why this factor is so large.) 

The maximum electric field is smaller for cores that collapse into larger black 
holes (because the discharge time during infall is greater) , so for a very conservative 
upper limit on the electric field, we shall assume that the black hole that forms has 
one solar mass. Of course, we expect that the minimum mass of a black hole that 
forms astrophysically has significantly more than one solar mass, so the correspond- 
ing maximum electric field would be weaker (by a logarithmic factor such that the 
maximum pair production rate would go essentially inversely with the square of the 
mass of the black hole). That is, we are assuming that one solar mass is a very 
conservative lower limit on the mass of any black hole that forms astrophysically in 
the present universe (as distinct from, say, primordial black holes that might have 
formed much smaller in the very early universe; our calculations will not apply to 
those, but they will have by now had a very long time also to discharge and so would 
also not be expected to have significant charge today). 

The maximum electric field is also smaller the slower that the core collapses 
(giving more time for discharge), so again to get a very conservative upper limit on 
the electric field, we shall assume that the core falls in as fast as is astrophysically 
possible, which is free fall with zero binding energy. We assume that it is not 
astrophysically possible to have the spherical outer boundary of a collapsing core 
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moving inward with a velocity so high that it would have come from an unbound 
configuration with nonzero inward velocity at radial infinity, or that at smaller radii 
it could have been accelerated inward at greater than gravitational accelerations. 

Thirdly, the maximum electric field is of course smaller the smaller the initial 
charge on the collapsing core. To get the maximum electric field possible, we shall 
assume that the initial charge-to-mass ratio of the collapsing core is unity, the largest 
value possible for which the electrostatic repulsive forces do not overwhelm the 
total gravitational attractive forces on the entire core and prevent the core from 
collapsing (since we are excluding the possibility that the core is shot in from far 
away or is otherwise pushed inward by nongravitational forces, which is not at all 
astrophysically plausible). 

Now of course if the charge-to-mass ratio of the core really were unity, the core 
would not fall until there is some discharge. But merely to get a conservative upper 
limit on the electric field as the core collapses, we shall for simplicity assume that 
the initial charge equals the mass but ignore the electrostatic repulsion, so that the 
core nevertheless falls in at the rate it would from purely gravitational free fall from 
infinity, with no reduction of the free fall rate by the electrostatic repulsion of the 
charge. The gravitational effects of the electric field outside the core would also 
reduce the free fall time, but we shall ignore this effect as well and simply take 
the external gravitational field to be given by the vacuum spherically symmetric 
Schwarzschild metric. 

First we shall discuss the physical quantities involved, quote the formula for 
the pair production rate, and cite the definition of a dyadosphere. Next, we derive 
a crude but very tiny upper limit on the maximum pair production rate of any 
macroscopic process and then make a rough estimate of the discharge rate for a 
collapsing charged core. Finally, we shall derive and solve the differential equations 
describing the process more precisely, in order to give the maximum electric field 
and pair production rate and the efficiency of the process. We also confirm that the 
interactions between the pairs produces are utterly negligible. 

2.1 Planck units and physical quantities 

Since the process involves both gravity (for the free-fall collapse of the core) and 
quantum field theory (for the discharge by Schwinger pair production), it is simplest 
to use Planck units, in which 
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G = k = C = 4:7l€o = /^Boltzmann = 1- (3) 

We can then always return to conventional units by using dimensional analysis 
to insert the right powers of the Planck units that we are setting equal to unity: 



1 = l p = JhG/c 3 « 1.61624 x 1(T 33 



cm, 



1 = t P = ^JhG/c 5 « 5.39120 x 1(T 44 s, 
1 = mp = Jhc/G « 2.17645 x 10~ 5 g, 



1 = S p = Jhc 5 /G « 1.22090 x 10 2S eV « 1.96561 x 10 ib ergs « 543 kwh, 



1 = T P = ^hc*/{Gkl olt _) « 1.41679 x 10 32 K, 

1 = F P = c 4 /G » 1.2103 x 10 49 dynes, 

1 = Np = Iphp 1 = c 7 /(h 2 G 2 ) « 4.39333 x 10 141 cm" 3 s" 1 , 

1 = p P = E P //p = c 7 /(hG 2 ) w 4.6331 x 10 114 ergs/cm 3 w 0.84960 x 10 200 eV/Mpc 3 



1 = Q P = v /47re /ic w 1.87554592 x 10" 18 Coulomb = 5.62274520 x 10~ 9 esu 
1 = R p = l/(47re c) = 29.9792458 ohms 
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1 = E P = ^c 7 /(A7re hG 2 ) w 6.4529 x 10 59 V/cm. (4) 

The numbers above were calculated from values given in the Particle Physics Booklet 
|3*Tj . When the ~ sign is used, usually the last two digits of the corresponding 
number are uncertain, but when the = sign is used, the result is exact (by definition), 
such as the Planck resistance Rp in ohms, which is 10~ 9 times the speed of light 
in cgs units, from the definition of the permeability of free space in mks units as 
fiQ = Att x 10 -7 newtons per square ampere. 

The results of our calculations depend mainly on just the following three physical 
quantities: the magnitude of the charge of the electron, positron, and proton, which 
we shall here denote by q (instead of the usual e, in order to avoid confusing it with 
the base of the natural logarithms); the mass of the electron or positron, m; and the 
mass of the sun, M & . For comparison we also give the mass of the proton, m p , and 
the mass of the neutron, m n . In Planck units, these various quantities are 

a w 0.08542454312, 
m w 4.18543 x 10~ 23 , 
M w 9.13617 x 10 37 , 
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m 



v 



7.68509 x 10" 20 , 

7.69568 x 10" 20 . (5) 



We shall find that the pair production process, for a maximally charged core 
collapsing to form a solar mass black hole, depends mainly upon the two large 
numbers I shall call A and B (and/or on a third rather large number C that depends 
logarithmically upon A and B and is independent of the solar mass M Q ): 

q 2 m 2 M 2 a ( GM G) m\ 2 , 
^ = - ° ss 3.39643251 x 10 2 

7T IX \ TlC J 



A = 



-.28 



5 = ? = J_f_?£_V/_^ ^ 42475 

" ~" 47rm 2 M 4tt \GM Q m) V 47re G'm 2 



ct / q 2 



.16.3 V4 , eoGm2 , -^.00843. (6) 

A is thus 7r times the fine structure constant a = g 2 / (47re ^c) ps 0.007297352568 
1/137.03599911 times the square of the ratio of the Schwarzschild radius of a solar 
mass black hole, 2M = 2GM Q /c 2 « 295 325.008 cm, to the Compton wavelength 
of an electron, A = 2?r/m = 2nh/(mc) w 2.426310238 x 10~ 10 cm. As we shall 
see, the largeness of A is the key to how much weaker than dyadosphere values are 
the maximum electric field and the maximum pair production rate that can develop 
from a collapsing charged core. (The experimental determination of A does not 
depend upon the relatively large uncertainty in Newton's constant G in cgs units, 
because it is GM & rather than M Q that is determined directly by Solar System 
measurements, so the value of A is known much more precisely than quantities like 
B that do involve the experimental determination of G.) 

If we use the well-known fact that in Planck units, stellar masses (particularly 
maximum white dwarf and neutron star masses) are of the order of the inverse 
square of the proton mass, 

M = m~ 2 ^ 1.69318 x 10 38 w 1.85327M© (7) 

(which coincidentally is within half of one percent of the largest prime found with- 
out computers, 2 127 - 1 = 170 141 183 460 469 231731687303 715 884105 727, which 
might be useful mnemonically to some readers), then one may give crude approxi- 
mations for A and for B that just depend on the elementary charge q and on the 
masses m and m p of the electron and proton respectively: 

~ q 2 m 2 M a ( m\ 2 ( fie \ a ( m\ 2 ~ „ _ _ , „™ 

A = z = - — -— = - — M 1.16654 x 10 29 , 

n 7i \m p J \Gm£J 7i \m p J 
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Thus A, and hence A, is very roughly a solar mass in Planck units, reduced by one 
power of the fine structure constant a and by two powers of the ratio m/m p of the 
electron mass to the proton mass. Similarly, B, and hence roughly B, is mainly the 
square of the ratio m p /m of the proton mass to the electron mass, but reduced by 
the square root of the fine structure constant and by the ratio of the area of a unit 
square to the area of a unit sphere. 

As one may derive from the definition below of a dyadosphere, an extreme 
charged Reissner-Nordstrom black hole (Q = M) has a dyadosphere if its mass 
is less than the critical mass 



so B is 1/4tt times the ratio of the critical mass M c to the solar mass M & . Later 
we shall see that it is impossible to form a black hole with Q/M ~ 0.99 if M < 
BCMq « 3 695 65OM ~ 3.7 x 1O 6 M . 

The combination C of the two independent large numbers A and B depends only 
on the properties of the electron and is the logarithm of the product of the reciprocal 
of 167r 3 , the fine structure constant a = q 2 = q 2 / (Aneohc), and the square of the 
magnitude q/m = (q/m)/y/47re G of the charge-to-mass ratio of the electron. C is 
large because the square of the magnitude of the charge-to-mass ratio of the electron, 
the ratio of the electrostatic repulsion to the gravitational attraction between two 
electrons, is very large: (q/m) 2 « 4.1657 x 10 42 . 

2.2 Pair production rate and definition of a dyadosphere 

The Schwinger pair production [29] gives a rate M of electron-positron pairs per four- 
volume (per three-volume and per time) in a uniform electromagnetic field that is 



where E and B are the electric and magnetic field magnitudes in a frame in which 
the electric and magnetic fields are parallel (which is a frame in which the electric 
field has the minimum value). 

Note that the exact one-loop pair production rate given above is not quite the 
same as twice the imaginary part of the one-loop effective action per four-volume, 



M c = q/m 2 = AnBM Q « 5.33751 x 10 5 M ( 



(9) 



[32j 




(10) 
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which is given by an infinite series J29J |o2] , of which the pair production rate is only 
the first term. However, for weak fields the first term dominates, so in such cases it 
is a good approximation to say that the pair production rate is twice the imaginary 
part of the one-loop effective action per four-volume. 

Here we shall be mainly interested in the case in which there is no magnetic field, 
in which case the pair production rate may be obtained from the formula above by 
taking the limit that B/E goes to zero: 

q 2 E 2 ( 7rm 2 \ m 4 e~ w 

47T 3 \ QE j 47T w 2 

where 

(12) 

if one defines the critical electric field strength E c to be 

2 2 3 

E c = — = « 2.05068 x 10" 42 « 1.32328539 x 10 16 V/cm. (13) 

q hq / \ J 

Ruffini [21 El lU 13 E] has defined a dyadosphere to be a region of spacetime in 
which the electric field E is greater than the critical electric field E c . Therefore, it 
is a region in which w < w c = n, and a region in which the electron-positron pair 
production rate is 

4 — 7T 

N > N c = t ~ 0.0003484287930 m 4 « 1.06924 x 10~ 93 (14) 

47T 7T^ 

« 4.69752273 x 10 48 cm" 3 s" 1 w 7.449533923 x 10 58 M" 4 . 

The last number means that in a cube of edge length M & ~ 1.47662504 kilo- 
meters, half the Schwarzschild radius of the sun, and in a time M & ~ 4.92549095 
microseconds, which would form a four- volume of M 4 w 1.58584308 x 10 10 cm 3 s that 
might be regarded as a minimum for typical macroscopic astrophysical four- volumes 
(about the same four- volume as a liter multiplied by half a year), more than 7 x 10 58 
pairs would be produced in a dyadosphere (under the highly hypothetical situation 
in which the dyadosphere could be maintained over this spacetime region, and also 
ignoring the effects of the pairs produced, as well as spacetime curvature). 

To put it another way, in a volume contained within a sphere in flat space of 
the Schwarzschild radius of the sun, 167rM| ) /3 ~ 53.9460048 cubic kilometers, the 
dyadosphere pair production rate would be more than 2.5 x 10 65 pairs produced per 
second. For the positrons produced to have a total charge equal to a solar mass 
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(in Planck units, as always when the units are not otherwise specified), one would 
need M Q /q w 1.06950 x 10 39 positrons produced, and the time it would take for this 
would be less than 10~ 26 seconds (again making the hypothetical assumption that 
the dyadosphere could be maintained over this volume for this length of time). 

This shows that the minimum dyadosphere pair production rate, though only 
47re _7r ~ 0.5430421124 in units of the electron Compton wavelength 
2vr/m Rj 2.42631024 x 10" 10 cm and the corresponding time of 8.093299792 x 10~ 21 
seconds, is utterly enormous at astrophysical scales. This strongly suggests that dya- 
dospheres will never form over macroscopic astrophysical scales. Indeed, the calcu- 
lations below confirm this and show that under extremely conservative assumptions, 
the pair production rate from astrophysical gravitational collapse is always less than 
10 -26 of that of a dyadosphere. (Under the plausible but not rigorous arguments of 
the Introduction that the charge is ejected when the electrostatic repulsion exceeds 
the gravitational attraction, the upper limit in a region of the order of the size of a 
black hole would be trillions of orders of magnitude weaker than a dyadosphere, a 
factor of more than 10 1Ql3 , and hence completely negligible.) 

2.3 Crude upper limit of pair production rate 

One can make the following crude estimate of the upper limit of the pair production 
rate from any macroscopic astrophysical process that has length scale L and time 
scale T. 

Suppose the electric field is 

p _ tiE c _ Tim 2 Q 

E = = = — , 15 

w qw L z 

thus determining a characteristic effective charge Q for the field producing the pairs 
at the rate per 4- volume given by Eq. fTTf . 



N = — N c « 1.07286 x 10 57 — cm^s- 1 . (16) 

w 2 w 2 

This gives 

w + 2\nw = \n(n 2 e 7T ^ > 1 (17) 

for M < ir 2 e 7T Af c = m 4 /(47r), which will be assumed for now and confirmed later by 
Eq. flUD below for T ~ L ~ 4.3 x 10~ 8 cm. Thus Eq. (JTHJ) gives Q < ixm 2 L 2 /q. 



11 



For the timescale of the electric field to be at least T, the charge Q' produced in 
this same time, which tends to discharge the electric field, cannot greatly exceed Q. 
This then gives 

tf~^ = £^-T*«<^, < 18 > 

or 

£ 5 J^ 8 . 45xl0 -§, (19 ) 

Al c q' L m z L± LI 
which is enormously smaller than unity for macroscopic astrophysical length and 
time scales L and T. 

This strongly argues that no matter what the situation, if the time and length 
scales obey 

27r 2 tt/2 

T ~ L > « 2.66 x 10 25 ss 4.3 x 10" 8 cm, (20) 

qm 

then dyadospheres of such a macroscopic size will not form. 

In the example below, we shall find a more precise upper bound on the pair 
production rate in the case of spherical symmetry that is slightly more than one order 
of magnitude stronger. This example confirms the limit above and strengthens the 
evidence that pair production rates in macroscopic processes must be many, many 
orders of magnitude below dyadosphere rates. 

2.4 Approximate estimate of the discharge rate 

We consider a positively charged spherical core of mass M freely collapsing into 
a black hole, with the surface at a radius R(t), and with vacuum outside (except 
for the electromagnetic field and pairs produced by it, which we assume will not 
significantly modify the Schwarzschild geometry; including such modifications would 
reduce the electric field even further than the conservative upper limits we shall find 
below). Outside the surface of the core, pairs will be produced, with the positrons 
moving outward and the electrons moving inward. (Since the core is expected to 
be highly conducting, virtually all of its excess charge would be at the surface, so 
inside the core there would be a negligible macroscopic field and hence negligible pair 
production.) We shall confirm in Section 4 that the interactions between individual 
positrons and electrons (e.g., the annihilation probability for each particle) is utterly 
negligible. 
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As the electrons pass inward through the outer boundary of the core over time, 
they reduce the value of the charge of the core, Q(t), limiting the value of the electric 
field outside, E(t,r) ~ Q(t)/r 2 , under the assumption (which can be verified from 
the results) that at any one time the charge contributed by the pairs outside the core 
radius R(t) is a small fraction of Q(t), so long as one does not go to such huge radii 
that it includes the outgoing positrons emitted over a large fraction of the previous 
infall of the core. (At such huge radii where the charge may differ significantly from 
that of the core at the same time, because of the outgoing positrons in between, the 
inverse-square falloff of the electric field would make it so weak that pair production 
would be negligible, so we need not consider such regions.) 

We also make the assumptions, which will be verified later (at least for M <C 
10 6 M Q ) that Q(t) <C M once the core gets within a few orders of magnitude of the 
Schwarzschild radius 2M, that the total energy that goes into the pairs is also much 
smaller than M (so that the core mass M stays very nearly constant), and that 
(with no restriction on M) the probability for any one of the particles produced to 
annihilate with an antiparticle is negligible. 

Because the pair production rate per four- volume, M = (m 4 / '4ii)w~ 2 e~ w , de- 
creases exponentially rapidly with w(t,r) = irm 2 /qE m 7rm 2 r 2 /qQ(t), and because 
the pairs produced decrease Q(t) and hence increase w(t,r) at fixed r, there will be 
a self-regulation of w(t,r) (described more precisely by the differential equations of 
the following subsections, though we do not need this for the approximate estimate 
of this subsection). In particular, if we define 

, 7iE c 7im 2 R(t) 2 . . 

which gives the pair production rate Af(t) = (m A /A-n)z~ 2 e~ z at the surface of the 
core (where the rate is maximal at that time, since the electric field has an inverse-r 2 
falloff at greater radii r > R), the self-regulation will keep z(t) changing only slowly 
with t and R(t). Hence Q(t) will vary roughly in proportion to R(i) 2 . 

This means that the logarithmic rate of change of Q(t) will be roughly twice the 
logarithmic rate of change of R(t). If we do a Newtonian analysis of free fall from 
rest at infinity, we find that dR/dt = —^2M/R, which makes the logarithmic rate 
of change of R(t) equal to — 1/(2M) times (2M/R) 3 ^ 2 . As the core approaches the 
horizon at R = 2M, this Newtonian estimate of the logarithmic rate will approach 
— 1/(2M), so the logarithmic rate of change of Q(t) will approximately approach 
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— 1/M. (Here we are ignoring relativistic corrections, but it turns out that they 
make only a relatively small difference in the maximum electric field.) 

Now we can calculate dQ/dt from the pair production rate (as a function of Q) 
and set it equal to —Q/M when R ~ 2M to solve for Q and hence the pair production 
rate. Because the logarithmic rate of change of Q is enormously less than that given 
above for a dyadosphere, the pair production rate must be suppressed by rather 
large values of w outside the core. Since w ~ -nm 2 r 2 j qQ increases proportional 
to the square of the radius, the suppression will rather rapidly increase with the 
radius outside the core (which we are now taking to be at R ~ 2M). This means 
that the pair production rate will decrease roughly exponentially with the radius 
and will reach a value a factor of 1/e smaller than at the core surface itself (where 
w = z = irm 2 R 2 /qQ ~ 47rm 2 M 2 /qQ) roughly when w increases by unity. Since the 
logarithmic rate of increase of w with r at r = R is 2/R, the point at which the pair 
production rate will have dropped by the factor of 1/e will be at r—R m R/ (2z) <C R, 
the last inequality coming because z ^> 1 in order that the pair production rate be 
much less than dyadosphere rates. 

With a roughly exponential decrease in the pair production rate with radius, 
at a logarithmic rate greater by a factor of z/2 ^> 1 than the logarithmic rate at 
which the radius grows, the total pair production rate per time is roughly the pair 
production rate per four- volume at the surface of the core (J\f(R) = (m A /4iT)z~ 2 e~ z ), 
multiplied by the effective volume where most of the pair production is occurring, 
which in this case is roughly the area 4nR 2 of the core surface (or of the black hole 
formed by the collapsing core) multiplied by the radial distance r — R « Rj (2z) 
out to where the pair production rate per four-volume has decreased by a factor of 
1/e. (For the pair production rate per coordinate time t, the relativistic correction 
that makes the proper radial distance greater than r — R is compensated in the 
approximately Schwarzschild metric outside the core by the relativistic correction 
that makes the proper time smaller than At by what is precisely the same factor in 
the Schwarzschild metric.) 

When the number rate per time is multiplied by the charge — q of each ingoing 
electron, one gets at R ~ 2M that 

M " -^W^Tz ~ ~P?— (22) 
For the logarithmic rate of change of Q to be roughly —1/M, we set dQ/dt equal to 
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-Q/M = —Trm 2 R 2 /(qMz) ~ — 4irm 2 M/(qz). This leads to the equation for z, 



■ 2 ~ n 2 M 2 ' ' r 1 

~K ~~ " I Mr. 



Z V ^ Z = q mM =M^\ = V ~ 3.39643251 x 10 28 /i 2 , (23) 



where 

„ = >L (24) 

is the ratio of the mass M of the freely collapsing core to the mass M & of the sun. 
Since A is so large, for |21n/i| <C In A one gets the crude explicit solution 

z ~ In A - 2 In In A + 2 In // « 57.325 + 2 In //. (25) 

One can also solve Eq. ()23|) numerically for z when \i = 1 (M = M©) to get what I 
shall call z*: 

z* « 57.588464. (26) 
Then for any |2 ln/x| <C In A, one gets the approximate solution 

z^z* + ^^-ln^ « 57.588464+ 1.932873 In /i. (27) 
z* ~\~ 2 

This means that the ratio of the electric field value E(2M) at the surface of the 
collapsing core (when it enters the black hole) to the critical electric field value E c 
for the definition of a dyadosphere is 

^1 = * « * - w 0.054552465 - 0.001830974 In/,. (28) 

Thus this approximate estimate would indicate that the electric field outside a 
charged collapsing core is always less than about 5.5% of that of a dyadosphere, 
differing from that of a dyadosphere by a factor of at least 18. (We shall see later 
that the numerical solution of the ordinary differential equation governing the charge 
during the collapse, including relativistic effects, agrees with the result above to 
about three decimal places, so the result above is quite accurate.) 

From this solution, one can see that the pair production rate at the surface of 
the core when it crosses the horizon at R = 2M is 

M { R) ~ £*V- ~^2 = 0^== i ^~7x 10-W. (29) 

which for M > M Q or \i = Mj M > 1 is more than 26 orders of magnitude smaller 
than the minimum pair production rate M c of a dyadosphere. 
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Therefore, the approximate algebraic estimate is that the pair production rate is 
never more than about 10 -26 times the minimum amount for a dyadosphere. In the 
next subsections we shall indeed confirm from the solutions of differential equations 
that this algebraic estimate is indeed a good approximation for the maximum electric 
field. This result implies that it is very unlikely that dyadospheres can form from the 
collapse of charged cores, even if somehow all discharge mechanisms are eliminated 
other than the pair production itself. 

2.5 Differential equations for the discharge 

In this subsection we shall derive and in the next two subsections solve (approx- 
imately) the differential equations for the discharge of the collapsing core. The 
differential equations will be derived under the assumption that the tunneling dis- 
tance for a pair to come into real existence, L tnnnc \ ~ mj (qE) = z/(nm) ~ 7 x 10~ 10 
cm, is much less than the astrophysical length scales for the collapsing core, which 
is a very good approximation. The differential equations will be solved under the 
approximation that this tunneling length is also significantly greater than the Comp- 
ton wavelength of an electron, which implies that w = TTm 2 /(qE) ^> 1 and hence 
that the pair production is mostly confined to a radial region r — R ~ R/w <C R 
that is much smaller than the radius R of the surface of the charged collapsing core. 
This latter approximation is less good but still leads to errors of only a few percent. 

After the charged particles are produced in pairs, they will be accelerated by the 
radial electric field (electrons inward and positrons outward, under the assumption 
here that the core is positively charged). Each time a charged particle travels a 
distance m/(qE) parallel to the electric field, it will gain an additional kinetic energy 
equal to its rest mass. Thus it will very quickly accelerate to a huge gamma factor 
and move very nearly at the speed of light. As a result, one will get effectively a 
null 4- vector of positive current density (highly relativistic positrons) moving radially 
outward and a null 4- vector of negative current density (highly relativistic electrons) 
moving radially inward. 

It is most convenient to describe this current in terms of radial null coordinates, 
say U and V, so that the approximately Schwarzschild metric outside the collapsing 
core may be written as 

ds 2 = -e 2a dUdV + r 2 {U, V)(d6 2 + sin 2 ^0 2 ). (30) 

For example, if one defines the radial variables u and v (not to be confused with the 
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null coordinates U and V above) and r* by 



u = v 2 = 



2M 



dr 



(31) 

i 

r = / - 2M \- ■ In (- l)l , (32) 

1 — u Yu \u J . 

then a common choice of null coordinates for the exterior region of the Schwarzschild 
metric is 



U = t 





1 


+ ln 




= t-2M 


_u 


\u 


= t + r* 










1 


+ ln 




= t + 2M 


,u 





- 1 



(33) 

U /J 

With this choice of null coordinates, one gets e 2a = 1 — u = 1 — v 2 . Other choices of 
radial null coordinates would be any smooth monotonic function of U as the retarded 
radial null coordinate, and any smooth monotonic function of V as the advanced 
radial null coordinate. Below we shall give another choice of null coordinates (U, V) 
that is useful for getting an approximation to the discharge rate. 

Now we can write the nearly-null outward number flux 4-vector of positrons 
as n + = n^dy and the nearly-null inward number flux 4-vector of electrons as 
n_ = n^du- Since each positron has charge q and each electron has charge —q, the 
total current density 4-vector is 



J = 



qui- 



qn\dy — qnFdu- 



(34) 



The radial electric field of magnitude E 
tensor 



Q/r 2 has the electromagnetic field 



F = --Ee 2a dU A dV = --^e 2a dU A dV, 
2 2r 2 



(35) 



where Q = Q(U, V) is the charge inside the sphere labeled by (U, V) (and which 
is a function only of these two null coordinates, because of the assumed spherical 
symmetry). 

Then from Maxwell's equations (essentially just Gauss's law here), we may de- 
duce that the null components of the current density vector are 

-2e- 2 °Q :U _ +Q> V 



■v v 
J = qn + = 



—qn_ 



4nr 2 
+2e- 2 °Q. 
4nr 2 



v 



Atct 2 ' 
_ -Q> u 
4nr 2 



(36) 
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Although of course the current density 4-vector field is conserved, the number 
flux 4- vectors n + and n_ of the positrons and electrons are not. Their 4-divergences 
are each equal to the pair production rate J\f when we can neglect annihilations, as we 
can here with the density of pairs being sufficiently small. When these 4-divergences 
are written in terms of the charge Q(U, V), one gets the following partial differential 
equation for the pair production and discharge process: 

^■-VW^^f-^f). (37) 

In covariant notation, with 2 OQ = —4e~ 2r7 Qw being the covariant Laplacian 
in the 2-dimensional metric 2 ds 2 = —e 2a dUdV and with □ being the covariant 
Laplacian in the full 4-dimensional metric, the differential equations for the pair 
production and discharge may be written as 

= 2 UQ = UQ- -Vr ■ VQ 

r 

= T U (Qj - rQD = rU (rE) - r 3 EO Q . (38) 

or, explicitly in terms of the electric field E = Q/r 2 , as 

9l . 2q 3 r 2 E 2 ( itm 2 \ , s , /1\ /N 

8TTqr 2 Af = -^— 2 — exp ( — —J = ^(r 2 ^) = rU (rE) - r ED i^-j . (39) 

Since for the sake of argument in this section we assume that the positively 
charged particles at the surface of the collapsing core do not escape to the outside, 
and since there is no electric field inside to produce particles there, the boundary 
condition at the surface of the collapsing core is that there is no outward flux of 
positrons there, so that Q jj = at the core surface. The boundary condition at 
infinite radii is that we assume that there are no incoming electrons there, but only 
outgoing positrons pair-produced by the electric field at finite radius, so at radial 
infinity, Qy = 0. The boundary condition in the infinite past is that we assume 
that Q is as large as it can be and still have the core collapse gravitationally, so 
there we set Q = M but then make the idealized assumption of ignoring the effect 
of the electrostatic repulsion on the collapse of the core. 

For an analysis of the solution of the partial differential equation (J37j) or (J38)) 
with the boundary conditions given above and with the surface of the collapsing core 
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freely falling inward from infinity in an assumed external Schwarzschild geometry 
(ignoring the gravitational effect of the electric field, which will be negligible once 
the core at large radius discharges enough to get its charge Q to become much less 
than its mass M), it is useful to shift from the commonly used null coordinates 
(U, V), given by Eq. ffTTTTf) above, to new null coordinates (U, V) which on the core 
surface are both set equal to the infall velocity of the core surface (as seen in an 
orthonormal frame carried by a static observer), v = \J2M/r, which runs from 0, 
at the idealized beginning of the core surface collapse at radial infinity, to 1 when 
the core surface crosses the event horizon at the Schwarzschild radius r = 2M and 
enters the black hole. (The variable v , and hence U and V at the core surface, would 
increase to infinity at the black hole singularity at r = 0, but since nothing inside 
the black hole can be visible to the outside, we shall ignore what happens inside and 
just consider the values of the electric field outside, where v < 1.) 

One can readily calculate that the freely collapsing core surface, moving along 
radial geodesies in the Schwarzschild metric with unit conserved energy-at-infinity 
per rest mass, has Schwarzschild coordinates 

t = 2M[- - V 3 - 2v~ l - In (l-t») + ln(l +v)}, 
2M 

r = — 40 

where we have chosen the origin of the t coordinate to avoid the explicit appearance 
of a constant of integration. The tortoise radial coordinate r* and the proper time 
r along the worldlines of the core surface have the form 

r* = 2M[tT 2 -21nt; + ln(l-i;) + ln(l + t;)], 

4M , . 

T = (41) 

where we have chosen the origin of the proper time to be the proper time at which 
the core surface hits the singularity at r = inside the black hole (so that r is then 
negative for all the physical times before then). 

Therefore, at the core surface at each value of v, Eq. fjHHjl gives the values of the 
null coordinates above as 

U = 4M[-- V - 3 --v~ 2 -v" 1 + \nv -\n(l-v)], 
3 2 

V = 4M[-i*r 3 + \v~ 2 - v- 1 -lnt> + ln(l +«)]. (42) 
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Now we want new null coordinates U = U(U) and V = V(V) both to equal v 
along the collapsing core surface, so the inverse relations must be 

U(U) = 4M[--U~ 3 - \u~ 2 -U' 1 +YnU -hx{l-U)\, 
3 2 

V(V) = AM[--V- 3 + -V- 2 -V- 1 -\nV + In (1 + 1/)]. (43) 
3 2 

The differential forms of these relations are 

AMdU 



dU 



u\i-uy 



AMdV 

dv = v^Tvy (44) 

Unfortunately, it seems unlikely that the direct relations U = U (U) and V = V(V) 
can be written explicitly in closed form in terms of elementary functions. 

In terms of these null coordinates (U,V), the Schwarzschild metric Eq. ffHOj) 
becomes 

, , 16M 2 (1 - v 2 )dUdV 4M 2 r ,„ 2 , „ , ,, 

ds - - uHi - u)vi + v) + — {de + Sln m ) ' (45) 

where v(U, V) is defined implicitly by 

v -i-2\nv + \n(l-v 2 ) = — = ^--^ 
V -2.\l 1.1/ 

= -u- 3 - -V~ 3 + -U~ 2 + -V~ 2 + u- 1 - v- 1 
3 3 2 2 

-\nU -lnV + ln(l-U)+ln(l + V). (46) 
It is convenient to define new time and radial coordinates 

T = 

X = (47) 



Then on the collapsing core surface, T = v = y2M/r = [4M/(— 3T)] 1 / 3 and X = 0, 
with — r being the proper time remaining until the surface reaches the curvature 
singularity at r = (where T increases to infinity, though T increases only from 
0, at t = — oo where the idealized freely falling core surface has r = oo and zero 
velocity, to 1 before the collapsing core surface enters the black hole horizon at the 
Schwarzschild radius r = 2M and the inward radial velocity v approaches the speed 
of light with respect to a sequence of static observers outside the black hole). X 
rises above as one moves radially outward from the collapsing core surface. 
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Nearly all of the pair production will occur close to the collapsing core sur- 
face (since the inverse- square falloff of the electric field with r will make the pair- 
production rate drop rapidly with r), where X < T 2 < T. Thus we can express the 
metric functions exactly as functions of T but for X <C T 2 < T only as an expansion 
in X/T 2 and in X/T. To second order in X/T 2 and in X/T, Eq. (|47)|) leads to 



T 



2T 3 



5T 2 

bl X 2 



u + v 



l^ + ^m^ (V-Uf. (48) 



Inserting this back into the metric Eq. (JSJ), and also substituting U = T — X and 
V = T + X, then gives the Schwarzschild metric near the freely collapsing core 
surface, to second order in X/T, as 



ds 2 



+ 



16M 2 

^8 

AM 2 



1 + 



4X2 

J^2 



-dT 2 + dX z 



AX 4 + 10T 2 \ 

X 2 (d6 2 + sm 2 



(49) 



For a good estimate of the dominant part of the pair-production rate and of 
the maximum value attained by the electric field outside the black hole, it seems to 
suffice to use this metric Eq. (J49j) simply to first order in X (which, when small, 
is proportional to the distance from the collapsing core surface, though with a In- 
dependent proportionality factor, so that the physical distance, in a radial spacetime 
direction orthogonal to the radially infalling geodesic worldlines along the collapsing 
core surface, is roughly AM/T 4 times X for sufficiently small X/T). If we return to 
the null coordinates (U, V) but keep only the single first-order term in IX = V — U, 
we get 



ds' 



2 V2 M 2 



dUdV + 



2 6 M 2 



(u + vy f 1 (u + v) 2 \ 



i + 



(d9 2 + sin 2 



(50) 



2.6 Approximate solution for the discharge 



In this subsection we shall construct an approximate solution of the partial differ- 
ential equation (|37j) for the pair production and discharge process. 

To simplify the constants appearing in the various equations, it is useful for each 
core mass M = M & n to define the normalized charge 



y 



q 3 Q 



An 2 Z Anm 2 M 2 ' 



qQ 



(51) 
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where q and m are the charge and mass of the positron (in Planck units) and where 
Z = q 2 m 2 M 2 /rr = Afi 2 was defined in Eq. flU, with A = q 2 m 2 Ml/n « 3.4 x 10 28 
being defined in Eq. (jUJ) and \x being defined in Eq. (J2l"j) as the ratio of the core 
mass M to the solar mass M & . 

It is also useful to use the quantity w = nE c /E = rrm 2 /(qE) = TTm 2 r 2 / (qQ) 
defined in Eq. (IT^j) and the quantity z defined in Eq. 12 which is the value of w 
on the core surface, where r = R. Then we have the simple relation 



between the radial variable v = y2M/r, the normalized inverse electric field w, and 
the normalized charge y. On the core surface, this becomes 



Of course, for any particular core collapse, M and hence Z are constants (since 
we are ignoring the tiny fraction of the total energy being carried away by the 
outgoing positrons and will later confirm that this fraction is indeed small, for 
M C 3 x 1O 6 M ). On the other hand, Q, E, Af, v, w, and y are scalar fields 
over the spacetime, functions of the spacetime point, depending on both the radial 
and time coordinates (t,r) or (T, X), or on the radial null coordinates (U,V) or 
(U,V), under the assumed spherical symmetry so that they do not depend on the 
angular coordinates (9,<p). 

An intermediate example is the normalized surface inverse electric field variable 
z, which is directly defined only on the core surface and depends only on the single 
variable time variable T — U — V — von that surface where this triple equality 
applies. However, below we shall often extend the definition of z to the region 
outside the core as either z(U) = w(U,U), the value of w at the core surface (where 
U — V) at the same value of U as the field point outside (where U < V), or as 
z(V) = w(V, V), the value of w at the core surface at the same value of V as the 
field point outside. That is, for a field point outside the core surface with radial 
null coordinates (U, V) (suppressing the irrelevant angular coordinates, since with 
the assumed spherical symmetry, all the fields we are interested in, such as the 
normalized inverse electric field w, are independent of them), z(U) is the value of 
w at the point on the core surface that is to the past of the field point, along the 
radial null geodesic that in the future direction of time goes radially outward, with 
constant U, from the point on the core surface where z(U) = w(U, U) is evaluated 




(52) 




T 4 zy = U 4 zy = V 4 zy = 1. 



(53) 
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to the field point at (U,V). Similarly, z(V) is the value of the normalized inverse 
electric field w at the point on the core surface that is to the future of the field point, 
along the other radial null geodesic that in the future direction of time goes radially 
inward, with constant V, from the field point (U, V) to the point on the core surface 
where z(V) = w(V, V) is evaluated. If the time coordinate T = (U + V)/2 is taken 
to run vertically upward, and the radial coordinate X — (V — U)/2 is taken to run 
horizontally to the right, so that the null coordinate U = T — X increases upward to 
the left and the other null coordinate V — T+X increases upward to the right, then 
the core surface is at the vertical line X = 0, and the three points (U,U), (U,V), 
and (V, V) form a right triangle with the right angle at the field point (U, V) outside 
the core, and the hypotenuse running vertically upward along the core surface (i.e., 
in the future direction of time along the surface) from (U,U) to (V, V). (Here of 
course U and V are the values they have at the field point, with U < V.) 

Using the normalized spacetime-dependent charge y = v ~ 4 w _1 instead of the 
charge in Planck units, Q = (Aimi 2 M 2 / q)y , the partial differential equation ()37|) for 
the pair production becomes 

42 -4-1 ze e 
V,uv = ~^Ve ' = (54) 

With the null coordinates U = T — X and V = T + X that were chosen to be equal 
to v on the collapsing core surface, inserting the exact and approximate values of 
e 2cr and of v into this differential equation gives 

y,uv = ^{y,TT-y,xx) 

8z(i - v 2 y y 2 



v y 



U A (l-U)V\l + V) 

%Z{l-v 2 ) e~ w 
U\l-U)V A {l + V)v 4 w 2 
128Z 



2 15 Z e~ l 
( U + v y2—2 

8Z e~ w 



^ TV , , (55) 
As with Q, y = qQ / \Airm 2 M 2 ) satisfies the boundary conditions yjj = at the 
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core surface (U = V), so that there is no outward flux of positrons there, and yy = 
at past null infinity (U = — oo), so that there is no inward flux of electrons there. 

To calculate the time rate of change of the normalized charge y at the surface, 
V,t = y,v there where y t u = 0, we can integrate y^jjydU radially inward along the 
inward null geodesic V = const, from U — at past null infinity to U — V at 
the core surface where we want to evaluate yy. However, to do this by Eq. 
we need to know y(U, V) all the way along this null geodesic, which is determined 
exactly by solving the nonlinear Eq. (jo3j) . with its boundary conditions, in the past 
wedge between the core surface and this inward null geodesic. 

Here I shall not attempt to determine the exact solution (which almost certainly 
would not be expressible in closed form in terms of elementary functions) or even 
to solve the partial differential equation (|55j) numerically (which will be left for 
an enterprising reader; I'm content to get a citation rather than another paper for 
this). Instead, I shall construct an approximate solution, first reducing the partial 
differential equation (Jo7)j) to an approximate ordinary differential equation for y(T) = 
y(v) or for z(T) = z(v) on the core surface, and then finding explicit approximate 
solutions and numerical solutions of this approximate ordinary differential equation. 

I shall use the fact that most of the pair production is confined to a region close 
to the surface, X < T 2 , where r = 2M/v 2 w (2M/T 2 )(1 + 2X/T 2 ) has not risen far 
above its surface value 2M/T 2 , so that the electric field has not dropped far below 
its surface value and so that the normalized inverse electric field variable w has not 
risen relatively far above its value z at the surface. (From the dominant factor e~ w 
in the pair production rate, the pair production becomes relatively negligible once 
w increases by a few units above its value z at the surface, but since z 1, w 
increases by a few units when its change is still relatively small in comparison with 
w or z.) 

Outside the dominant pair-production region, there may be a significant out- 
ward current flux qn + = qri^dy oc —yjjdy of positrons flowing outward from the 
pair-production region, but the inward current flux — qn_ = —qn^djj oc +yydjj will 
become negligible once one gets sufficiently far outside the dominant pair-production 
region (since there is assumed to be no current coming inward from infinity). There- 
fore, outside the dominant pair-production region, y will be nearly constant along 
the outward null geodesies, U = const.. Within the pair-production region, the 
inward number flux 4-vector n^du of electrons will be future directed, so nP will be 
positive and hence yy will be negative. This means that as one moves outward along 
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a outward null geodesic U = const, from the core surface, y will initially decrease 
a bit and then rapidly settle down to near its asymptotic value, the charge as eval- 
uated a long way away at fixed retarded time U. (For this one should avoid going 
so far away that the tiny difference from the speed of light of the velocities of the 
outgoing positrons becomes significant, but that would be at an enormously large 
radius. In this paper, unless explicitly stated otherwise, I am assuming that the 
electrons and positrons propagate effectively at the speed of light after their enor- 
mous electrostatic accelerations that give each of them an increase of their energy 
equal to their rest mass m each time they travel a radial distance w/(irm) <C r.) 

Therefore, a zeroth-order approximation would be to take y to be constant along 
the outgoing null lines, so that 

y(U, V) « Vl (U, V) = y(U, U) = lJ- 4 w(U, U)' 1 = U^z{U)-\ (56) 

Henceforth I shall assume that z = z(U) = w(U,U), the value of ttE c /E ee 
Tim 2 / {qE) at the core surface at the same retarded null coordinate U as at the field 
point (U, V) that is along an outward future-directed null geodesic from the point 
(U, U) on the core surface. Since the electric field decreases as one goes outward 
from the surface, whereas on the surface it increases with time as the core contracts, 
z is inversely proportional to the highest value of the electric field to the causal past 
of the field point (U, V) and is the lowest value of w to the causal past of (U, V). 

Now write the normalized charge y(U, V) as the zeroth-order part yi(U, V) that 
actually depends only on U, plus a remainder y 2 (U,V) that vanishes at the core 
surface U = V and which will generally be smaller than yi(U, V) by a factor of 
roughly 1/z: 

y(U, V) = Vl (U, V) + y 2 (U, V) = z^U^ + y 2 (U, V). (57) 

In the region near the core surface where the dominant pair production occurs, 
one can use Eqs. (j^Hj) and (jBTj) to write 

W (U, V) = v- 4 y- 1 ^z + Az^-^X - z 2 U 4 y 2 . (58) 

U z 

One then inserts this into Eq. (jB3j) and changes coordinates from (U, V) to (U, X) 
to obtain, for X = (1/2) (V -U) < [/, 



1 I 8Z e~ w 8Z 

' yi2 

8Ze 



^y2,XX + ^V2,XU ~ ji2 ^2 ~ U 12 Z 26 



U 12 Z< 



e -^X + ^U*y 2 _ (5Q) 
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Here I have used the fact that since y\ depends only on U, it gives yi t uv — 0, and 
the fact that the relatively small y 2 term makes a significant contribution on the 
right hand side only to the e~ w term. The effect of higher-order terms in y 2 in 
the exponential are negligible so long as z 3 U s y 2 <C 1, which turns out to be the 
case, except possibly very near the black hole horizon during the final stages of the 
externally visible part of the collapse. 

The boundary conditions of this partial differential equation for y 2 are that it 
vanishes along the core surface at X = 0, and asymptotically for large X one has 
yi,x = to reflect the fact that there is no incoming flux of electrons at large radius, 
outside the pair production region. 

Because z ^> 1 implies that the dominant pair production region occurs near the 
core surface, it turns out that the yi,xv term is generally negligible in comparison 
with the y 2 ,xx term (again except possibly very near the horizon of the black hole 
that forms as the end result of the gravitational collapse). Therefore, for each 
value of U and z = z(U), the partial differential equation (|59p becomes an ordinary 
differential equation for y 2 (X). Although the y 2 term in the exponent makes it 
a nonlinear ordinary differential equation, it is readily solvable when one chooses 
the exponent to be the dependent variable, with the solution (when added to y\ 
to get the total normalized charge y = qQ / (47rm 2 M 2 ) inside the sphere at radial 
coordinate X and retarded time coordinate U) being 



zU A I z 



-(v / r+p + i)--(v / r+p-i)e 



1, n. ~4z^fX 

u 2 



(60) 



Here 



= 4Ze- = 4C/(1 - S) 

~ u'(i - u) 2 z 2 ~ (i -uy ' y ' 

Since z = z(U) is independent of X, both P and S are functions only of the retarded 



where 



null coordinate U that was chosen to be the value of the infall speed v = ^2M/R of 
the core surface, at r = R(U), relative to an observer at fixed Schwarzschild radial 
coordinate r. 

In order that there be no outgoing flux of positrons at the surface (since there is 
no electric field inside the conducting core and hence no pair production there), one 
needs that the partial derivative of y with respect to U at fixed V be zero at the 
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core surface. In terms of the coordinates (U,X), this condition at X = becomes 

111 

= y,u- -jV,x = yi,u ~ ■zVi.x + V2,u - -jV2,x- (63) 

However, y\ has no dependence on X, so the second term vanishes, and along the 
surface at X = one has y<i = 0, so the third term vanishes as well. Therefore, this 
boundary condition relates the rate of change of the normalized charge evaluated 
at the surface, which is purely yi(U), to the radial derivative of yi at the surface 
(where y2 vanishes, but not its radial derivative). Then when one differentiates the 
relation z = U~ i y~ i along the surface, one gets 

- 8SZ/U (64) 



i + u + ^/(l + uy-Aus 

Because this equation is evaluated at the core surface (X = 0), the X-dependence 
has dropped out, and one thus has a nonlinear first-order ordinary differential equa- 
tion for z(U). The boundary condition for this equation is that at the infinite past 
(for the proper time r), when the core surface radius is at R = oo and hence has 



U = v = y2M/R = at the surface, we assumed the highest possible value for Q, 
namely Q = M, so then y(U) has the initial value y$ (at radial infinity) 

_ qQo q B 

m=m = Wi^ = WM 0/ , = 7 (65) 

where B w 42475 was defined in Eq. © and /x is M/M Q , the core mass in units of 
the solar mass M & . Thus initially z = U~ 4 y~ 1 ~ (fi/B)U~ 4 is infinite at U = v = 0, 
and S = 1 there, so that dz/dU ~ —Az/U is negative infinite initially. 

The Newtonian limit of Eq. (}64j) is the limit in which [/< 1, since U represents 



the inward velocity of the core surface that is at U = V = v = J2M/R, relative to 
an observer that stays at fixed r. In this limit, one gets 

dz . z ,z(* Ze~ z \ 

— « -4— S = -4— 1 — . (66) 

dU U U\ U 5 z 2 J y 1 

Since this is an ordinary differential equation along the worldline of the core surface, 



it is convenient to replace U by v = J 2M/ R, which is more obviously a coordinate 
along the core surface. Then Eq. ()66|) becomes 
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If one uses z = v ~ 4 y~ 1 , one can rewrite this explicit differential equation for z(v) as 
the following explicit differential equation for y(v) = yi{v), the normalized charge 
at the core surface: 

^K-AZvVe-"-**- 1 . (68) 
dv 

These differential equations, Eqs. (|fi7Jl and (jEBJ), are those of the discharge 
process in the Newtonian limit of the core collapse. They may be readily derived in 
that limit by assuming that when the electric field produces pairs of electrons and 
positrons, the electrons instantaneously propagate to the core to reduce its charge, 
and the positrons instantaneously propagate to large radii, so that the charge Q is 
at each moment of time independent of the radius. 

Actually, to derive Eqs. ()67|) and in the Newtonian limit, one assumes that 
z ^> 1, which leads to a simplified approximate expression for the integration of 
the pair production rate over all radii outside the core surface at each instant of 
time. Without this approximation, but still taking the Newtonian limit of the core 
surface moving very slowly compared to the speed of light, so that the electrons can 
propagate to the core and the positrons can propagate very far from it in a time in 
which the core radius changes only very little, one gets that the factors containing 
Z on the right hand sides should be multiplied by 

/DO / V \ 3 / 2 , 

e z ~ w - dw = 2z - 2v / ^e 2 erfc^ 

= 1 -^ + ^- i 8^ + ^ < 69 > 

a factor that is always just a bit smaller than unity (for the inevitably large z). 
Then Eqs. flBTj) and (J55J) become (after also multiplying the latter equation by v/y) 



,'^«-4(l-^W]=4(l-e--), (70) 

z dv \ v b z 2 



y' = -^^ -4Z V 3 3 /V(v~V 1 )e- w " V1 ee -4e" w , (71) 
y dv 



where 



W = z + 2\nz-\n J(z) -lnZ + 51nw. (72) 
Using the fact that J(z) obeys the differential equation 

'" 'l + (T3) 



dz V 2z 
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with the boundary condition that J(z) — > 1 for z — > oo, one can readily show that 
differentiating Eq. (J72|) gives 

so in this Newtonian analysis W cannot decrease to 0, and z simply continues to 
decrease as v increases. 

One may now insert the correction factor J(z) into the relativistic equations 
above to get improved equations that work to all orders in 1 j z when v is very small, 
though there will be error terms of order v/z that I shall not bother calculating. 
Again replacing U by its value v on the core surface where these ordinary differential 
equations apply, one can define 

S = l-^^ = l-e- w (75) 

and 



(™) 



AZe~ z J(z) = Av{l - S) 
v 4 (l - v) 2 z 2 = (1-v 

Then the relativistic versions of Eqs. (fTOj) and (fTTj) become 



v dz —8S 



zdv *i _i_ ii _i_ a fn^~ntf 



l + v + J{l + v) 2 -AvS 



(77) 



,_vdy , -8(1 - S) 

y = -— = -4- z w =■ (76 

1-^ + ^/(1 + ^)2-4^ 



(It is not a typographical error for the first denominator to have 1 + v outside the 
square root and for the other to have 1 — v; one can straightforwardly check that 
the two equations are consistent with each other.) 

The approximations used in my derivation of these equations are only valid for 
1 — v ^> 1/z, so they will break down as the core surface nears the black hole 
horizon, where v — 1. However, the ordinary differential equations (J77|) and (J78)) 
do not have any singular behavior as v approaches unity, so it might be that they 
remain fairly good approximations even as the core surface enters the black hole. 
On the other hand, I have definitely neglected various terms of the order of v/z, so 
near the horizon there will almost certainly be relative errors of the order of 1/z, 
which as we shall see below is of the order of 2%. 
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2.7 Approximate solutions of the ordinary differential equa- 
tions 



Now let us find explicit approximations for the solutions of these ordinary differential 
equations and compare the results to the numerical solutions. 

First, let us define various quantities that will be useful in the analysis. Some 
useful constants (depending on the mass M = M Q fx of the collapsing core) are 

(o 2 m 2 M 2 \ 
^— J « 65.69511 + 2 In //, (79) 

B 42475 

y = - « , 80 

/i // 

L = L + - lny = m A + - mL 5 + - In// « 79.01593 + 0.75 ln/x. (81) 
Second, it is convenient to remember the definition 

J{z) = 2z- 2v^IV erfcv^ (82) 

and also define some new explicit functions of z = nE c /E at the collapsing core 
surface, such as 

3 21 69 

K(z) = z + 2\nz -In J(z) = z + 2\nz + — - — + — + 0(^ 4 ), (83) 

its derivative with respect to z, 

rW , , dFsr(z) 1 2 3 21 207 5 , 

W^s) = tf(z) -^lnz-L, (85) 

and 

ee ee Z^'V^e"*./^). (86) 

Note that and are two different functions of z, though they are approximately 
the same in certain circumstances. 



Third, it is convenient to define several explicit functions of v = y2M/R (with 
R being the core radius during its assumed free-fall collapse from rest at radial 
infinity), such as 



/(u) = L -5 lm, = lii 1-^1, (87) 



30 



I \ — -1-4 

zi{v) = y v , 



41n/»- 0.25 + 1.25?; 4 In 2 /(v) - (8.5 - 2.5^) ln/» 
22(v) = /(u)-21n/(i;) + jt^ + 



/(«) f(v] 



2 



(89) 

X l (v) = exp [K(zl(v)) - f(v)} = l ) \ (90) 

5(1+^) 5(l + 2r;)(3-2r;) 
^ = ~4^(u) ' (91) 

X^X^+X^), (92) 
^) = /(„) + In [1 + X(t>)], (93) 

3 45 
z(v) = g{v)-2\ng{v) + g(v)-\A]ng(v) - -) + g(v)~ 2 (A\n 2 g(v) - ll\ng(v) + — ) 

+ g(v)^ In 3 «,(«) - 30 In 2 g(v) + y In <?(<;) - (94) 

^ s zw) ' (95) 

_ 47rm 2 M 2 y(t;) M o/ x 2 ^) 

= = 5 ' ( 96 ) 

q t> 

- , , irm 2 , , 

^ SM' (97) 

AA W - „ (98) 

and 

P «) s — , 99 

Finally, it is helpful to define or recall some functions that are explicitly given 
in terms of both v and z, though of course after solving the differential equation 
(177)1 for z(v) or for v(z), they can then be expressed as functions purely of v or 
of z, though there will not be closed-form explicit exact forms for that functional 
dependence: 

Y = e~ w = l-S= Ze ~l J y\ (100) 
W = -lnF = -ln(l - S) = K(z) - f(v), (101) 
X = ^ ee e w - 1 ee exp [A-(z) - /(«)] - 1 = ^ - 1, (102) 
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and 

2 



l + v + J(l-v) 2 + AvY 




(103) 



Here we shall motivate the choice of z(v) as an excellent approximation (for 
M < 3 x 10 6 M Q ) to the solution of the differential equation (J77j) for z(v), which 
then gives y(v) as an excellent approximation to the normalized charge y(v), Q(v) 
as an excellent approximation to the charge of the core surface when it has radius 
R = 2M/v 2 , E{y) as an excellent approximation to the electric field at the core 
surface, Af(v ) as an excellent approximation to the pair-production rate at the core 
surface, and p(v) as an excellent approximation to the ratio of the pair-production 
rate at the core surface to the minimum value in a hypothetical dyadosphere, which 
has E > E c = m 2 /q. 

From the definitions above and the differential equation (|77jl . one can readily 
derive that with the independent variable being chosen to be f(v) = L — 5\nv, 
X = exp [K(z) — f(v)] — 1 obeys the differential equation 

dX AK'zX , s 

T~ — ~ x (104) 

where here z is to be defined implicitly as a function of / and X by inverting the 
formula for X to get 

K(z) = / + ln(l + X) (105) 

and then solving this for z. For example, for K(z) ^> 1, which is always the case 
here, an approximate solution for z is 

z w K-2hiK + (AlnK--) + K~ 2 ( A In 2 K - 11 In K + 4) S 



+ K- 3 (— In 3 K - 30 In 2 K + — In K - — V (106) 
V 3 2 8 / 

Atf = 0or/ = L — 51nt> = oo, one has that AK'z = oo and D — 1, so X — oo 
there. One can see directly from Eq. (J78)) that the normalized charge y stays very 
near its asymptotic value y(0) = yo = B/fi for sufficiently small v (where the core 
is so large, and its electric field is so weak, that there is negligible pair production), 
so there z = y^v^ A w l/o" 1 ^ -4 = z\{v). Then X\(v) = exp [K(zl(v)) — f(v)] obeys 
the equation 

<JX ! ^-K'^zM - l) X lt (107) 



df V5 
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which is precisely the same as the Eq. (I1(J4|) that X obeys if D = 1, if z is replaced 
by z±(v), and if the final 1 is omitted from Eq. p()4jl . 

D«l when either X ^> 1 or v <C 1 (or both), so that one may then expand Eq. 
(1103)) for D to second-order in v to get 

v v 2 X 

D^l + ——+ , (108) 
1 + X (1 + X) 2 K 1 



Therefore, until v = y2M/R gets to be near unity (its value when the core surface 
collapses into the event horizon at R = 2M), it is a fairly good approximation to 
take D pa 1. 

The relation between Zi(v) and Xi(v) is 

K(z 1 (v)) = f(v)+\nX 1 (v), (109) 



which differs from Eq. (jl05j) in not having the additive term 1 in the argument of 
the logarithm, so X ~ Xi implies that \z — Z\\ <C 1 only if X\ ^> 1. However, until 
X\ drops to become much smaller than unity (e.g., for all X\(v) ^> ziiy)^ 1 <C 1), 
one still has z\ ^> 1, so as far as the relative error is concerned, one still has z pa z\ 
in Eq. (ITUl) . 

Finally, the final 1 in Eq. (|104|) is negligible so long as z\{y)X\{y) ^> 1, which 
for \i = M/M Q < B w 42475 is sufficient to imply that both D pa 1 and z\ > 1. 
Therefore, under these conditions (sufficiently small t>), -X"i(i>) is a good approxima- 
tion for X{y). Then one may use Eqs. ()105j) and (jl06|) to get a good approximation 
for z(v ) for these small values of v. 

Now let us discuss what happens when we go to larger values of v , so that zX 
becomes no longer large and Xi(v) is no longer a good approximation to X. If z 
and D are taken to be known functions of v (which of course is not really true before 
one finds the solution, but one can use crude approximations for them), then Eq. 
()104j) is a linear equation in X and hence has, as v increases and / decreases, X 
exponentially approaching a separatrix. Therefore, for v significantly larger than 
its value where zX becomes of order unity, X should have become very near the 
separatrix. Thus we can calculate an approximation for the separatrix in order to 
get an excellent estimate for what X (and hence z) is for large v (e.g., v = 1, where 
the core surface enters the black hole horizon). 

This behavior of the differential equation 1)104)1 is a consequence of the physical 
fact that the pair production is a highly regulated process. Once pair production 
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becomes significant to eject charge at a logarithmic rate comparable to that of the 
collapsing radius R(t) (which is far, far before dyadosphere values are reached, since 
the astrophysical length and time scales of the collapsing core are far greater than 
the Compton wavelength of the electrons and positrons being produced), the pair 
production keeps discharging the core as its radius contracts and prevents its electric 
field from rising very fast. 

As a first approximation, the right hand side of Eq. (jl(J4j) tends to zero, so that 
X tends to become nearly constant. But of course as z (and, to a lesser extent, 
D) changes, the right hand side of Eq. (jl()4)) cannot be precisely zero without X 
changing, so one really needs to keep both the left hand side and the right hand side 
slightly nonzero along the separatrix that X(f(v )) approaches. 

One can postulate that the separatrix can be written as a power series in z~ x 
(which stays small all the way to v — 1), with coefficients that are functions of v. 
One can then show that to second order in z~ l , the separatrix is at 

X(v, z) « -(1 + v)z- 1 - — (1 + 2u)(3 - 2v)z- 2 . (110) 
4 lo 

One can now insert this in place of X into Eq. (|l()5j) to obtain an equa- 
tion for z(f), after inverting the definition f(v) = L — 51nt> in Eq. (|87p to get 
v = exp [(L — /)/5] as a function of /. Because Eq. (jllOj) gives X(v,z) as an 
explicit function of both v (or of /) and of z, inserting it into Eq. (jl(J5|) will give 
^-dependence on both sides of that equation, so one cannot directly use Eq. ()106|) 
as an approximate solution for z(v) or z(f). However, one can find that an approx- 
imate solution for this large- v regime is given by the function z 2 [v ) defined by Eq. 
(JBHJ). If we then insert this approximation for z into Eq. (jllOj) for the approximate 
separatrix given X(v,z), then this becomes the explicit function X2(v) given by 
Eq. (|91|). This is then an excellent approximation to X(v) at large v, where the 
approximation valid at small v, Xi(v), becomes much smaller than 1/z. 

Thus we have produced approximations for both small and for large v , divided 
by a transition region where the right hand side of Eq. ()1U4|) is of the order of unity. 
However, we would ideally like a single formula that applies for all v , including the 
transition region. It turns out that a simple procedure works quite well, essentially 
because v , z, and D change very little across the transition region where the right 
hand side of Eq. ()104|) rapidly drops from being much larger than unity for smaller 
v to much smaller than unity for larger v. For smaller v, where Xi(v) is an ex- 
cellent approximation to X(v), Xi(v) ^> X 2 (v), but for larger v, where X 2 {v) is 
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an excellent approximation to X(v), ^(f) ^> X\(v). Therefore, in each case the 
better approximation dominates, and one can simply take their sum as the excellent 
approximation X(v) = Xi{v) + X 2 (v) given by Eq. 

When X(v) is inserted as an approximation to X in Eq. (|1U5|) . and when Eq. 
(jlOfij) is then used to get an explicit approximate solution of this equation for z(v), 
one obtains z(v) given by Eq. ()94j1 . where the g(v) used therein is defined by Eq. 
()93|) as what one would get for K(z) from Eq. f)105j) with X(v ) used in place of X. 
As a result, we have obtained a completely explicit approximation z(v) for z(v ) for 



all f = y2M/R from t> = (core surface radius R = oo) to t> = 1 (core surface 
radius R = 2M, crossing the black hole horizon). 

The only part of the explicit algorithm for getting this approximate expression 
z(v) that is not an elementary function is the definition of J(z) in Eq. (|82jl in terms 
of the complementary error function. However, an elementary function that gives 
an excellent approximation to J(z) for the large z that occur in the problem, with 
error only 0(l/z 5 ) that turns out to be numerically less than 2 x for all values 
of z that occur, is 



J{z) 



i + f~ + Ve1 z-i 



11 i 7 



y/6) 



z- 1 



11 _ 7 
'IS 75 



[111) 



If this formula is used instead of J(z) from Eq. (|H2*|) in Eq. (jHHj) for K(z) and 
their subsequent use in Eq. (jHUj) for defining X\(v) that then goes into Eq. for 
defining X(v), which then in turn is used in Eq. (}93|) to define g{v), which finally is 
inserted into Eq. (JMj) to get the final approximation z(v) for z(v), we then have a 
completely explicit elementary (though admittedly slightly involved) formula for an 
approximation to the solution to the relativistic ordinary differential equation (|77jl . 
which itself came from an approximate solution to the covariant partial differential 
equation (|3*H|) for the pair production and discharge throughout the spacetime region 
outside the surface of the collapsing core (assumed to be freely collapsing in with a 
Schwarzschild metric outside, ignoring the electrical and gravitational effects of the 
electric field, both which would slow down the collapse and increase the time for 
discharge, thereby increasing z(l) and reducing the pair production rate). 

If we would like a simpler approximate formula (though of course less accurate), 
we may drop the factors of J(z) (since for z > 57.6 it is always within 2.5% of 
unity), replace X 2 (v) by zero (since it is always small), and take only the first two 
terms of Eq. (|94jl . This then gives the approximation for z(v) that I shall call z(v), 
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z{v) = In ( V^" 5 + B~ V^~ 8 e B "^" 4 ) 

- 2 In In (A/i 2 ^ 5 + B~ 2 ^v^e 13 " 1 ^). (112) 

Once we have the approximate solution z(v) for z(v) = tcE c /E, where E is 
the electric field strength at the surface of the freely collapsing core when it has 
radius R = 2M/v 2 , we can easily evaluate other quantities on the core surface. For 
example, an approximation to the normalized value y of the charge Q is given by 
y(v) = l/[v 4 z(v)] from Eq. (J95)) . an approximation to the charge itself is given 
by Q(v) = 47tm 2 M 2 y(v)/q from Eq. (J96|) . an approximation to the electric field 
E is given by E(v) = rrm 2 /[qz(v)] from Eq. ()97|1. an approximation to the pair 
production rate per four- volume is given by J\ (v) = m 2 / (47rz 2 e z ) in Eq. ([98)1 . and an 
approximation to the ratio of this rate to the minimal dyadosphere rate M c (that of 
the critical electric field E c = m 2 /q) is given by p(v ) = Af(v)/Af c = exp (n — z)tt 2 /z 2 
in Eq. 

Of course, one can do the analogous procedure from the simpler approximate so- 
lution z{y) to define y{y) = l/[w 4 5(f)], Qiy) = 4Trm 2 M 2 y( y v)/q, E(v) = nm 2 /[qz(v )], 
J\f(v) = m 2 /(A7[z 2 e 2 ), and p(v) = N(v)/N c = exp (tt - z)tx 2 /z 2 . 

Now let us compare these explicit approximate expressions with the numerical 
solution to the ordinary differential equation (fTTj) . When this was solved numerically 
with Maple for 1 < fi < 3 x 10 6 (core mass M between one solar mass M & and three 
million solar masses), it was found that the more detailed approximate solution 
(barred quantities) agreed extremely well with the numerical solution over the entire 
range of v. For example, for ji — 1, it was found that the approximate solution ziv ) 
differed from the numerical solution z(v) by a relative error that was always less 
than 3.2 x 10 -5 ! The relative error had its maximum absolute value at the peak of a 
very small blip in the transition region at v ~ 0.0236144, where [z(v) — z(v)]/z(v) ~ 
-0.0000317. 

The relative error was very much smaller for v significantly lower than the transi- 
tion region, which is not surprising, since there the pair production rate is negligible, 
so that it is a very good approximation just to take the charge to have its initial 
value (which was Q = M at very large radii, but ignoring the effect of the energy 
density of the electric field on the metric and infall rate). For v above the transition 
region, the relative error was roughly constant and positive, but with a value sig- 
nificantly smaller than the magnitude of the small negative peak in the transition 
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region. A characteristic value of the relative error for the larger v region is given by 
the value at v — 1, when the core surface just enters the black hole horizon. 

A quadratic fit of the numerical solution for z(l) at the horizon (v = 1) versus 
fi for pL = 1 (M = M ), \i = y/lO, and /i = 10 (M = 1OM ) gave 

z(l) « 57.60480 + 1.932412 In /i + 0.0010297 In 2 /i, (113) 

whereas the explicit formula gave the quadratic fit 

z(l) « 57.60483 + 1.932414 In// + 0.0010281 In 2 (114) 

so the relative error in z at the horizon was about 5 x 10 -7 , or about one part in 
two million. Although this seems remarkably small, it is actually not too surprising, 
since X(f), and hence also z(v), rapidly approach the separatrix as v is increased 
beyond the transition region where the pair production starts to become important. 

One can see that the approximate estimate given in Eq. (|27jl in Section 2.3, before 
the differential equations for the pair production and core discharge were derived 
and solved, z(l) ~ 57.58846 + 1.932873 In /i, is also quite close to the numerical 
result, with a relative difference of only about —0.0003. Therefore, the approximate 
estimate made there was quite good, though its difference from the numerical result 
for z(l) is about 600 times the difference between the much more sophisticated 
estimate of z(l) and the numerical result for z(l). 

The simpler approximate solution z(v) of Eq. ()112j) has a somewhat larger error, 
though it always appears to be within 1% of the correct answer. The maximum 
absolute value of the relative error appears to occur at the black hole horizon (v = 
1) when fi « 3 693 360, where the numerical solution gives z(l) ~ 88.29447 but 
z(l) « 87.46451, with a relative error of about -0.940%. For /i < 3 539 000, the 
relative error of z(v) always seems to have magnitude always less than 0.5%. For 
example, for fi — 1 (M = M Q ), the maximum absolute value of the relative error is 
again at the horizon, v = 1, as it always seems to be for any fi, and at this value of 
(jl, the relative error of z{v) is about -0.486%. 

One might think that the simpler approximate solution is thus adequate, but 
because z is large and is exponentiated in the formula for the upper limit on the 
pair production rate Af, even a small relative error in z(v) can lead to a large relative 
error in J\f, as we shall discuss in more detail below after giving the numerical results 
for M . The result is that although z(v) is adequate for getting an estimate, within 
1% of the correct value, of z(v), and hence also of the upper limit of the electric 
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field E = TTm 2 /(qz) at the collapsing surface of the core of mass M = M®fi when 
it has radius R = 2M/v 2 , z(v) and the corresponding J\(v) are not adequate, even 
within 30% accuracy, for giving the upper limit on the pair production rate M . 

One can now also use the numerical results for fi = 1, fi = vTO, and /i = 10 
to calculate that for a solar mass core collapsing freely inward from rest at radial 
infinity with initially maximal charge Q = M (but assuming the Schwarzschild 
metric and infall rate), the (extremely conservative) maximal ratio of the electric 
field to that the critical field of a hypothetical dyadosphere, the value at the horizon 
at R = 2M, is 

E^Ml = _?L ~ 0.05453699 - 0.001824768 In u + 0.0000540770 In 2 fi. (115) 
E c z(l) 

The ratio corresponds to an upper limit on the electric field itself of 

-Emax = E(2M) (116) 

« 7.216801 x 10 14 V/cm 



( M \ , , ( M 

1-0.03345927 In — +0.000991565 In 2 



E(2M)/E C is again close to the estimate of Eq. (JH, 0.054552465-0.001830974 In fi. 
Thus we have confirmed that indeed the value of the electric field of a collapsing 
core is always less than 5.5% of that of a dyadosphere, differing from that of a 
dyadosphere by a factor that is always more than 18. Of course, in reality it would 
be unlikely ever to get a collapsing core with Q anywhere near M, so this is only an 
extremely conservative upper limit on the electric field. 

Now we can also use the numerical result to calculate an upper limit on the pair 
production rate Af, say as a ratio of that rate with the minimum pair production 
rate Af c of a dyadosphere, with fi = M/M Q : 



J\f(M) ttV 

Pi 1 ) = 



AT C z(l)V« 
6.61168 x 10~ 27 



/' 2 



1 + 0.0005545 In// - 0.00001759 In 2 fi). (117) 



Therefore, assuming that M = M & or fi = 1 is a very conservative lower limit on 
the mass of a core that can collapse into a black hole, we indeed see that even if one 
can somehow start with Q = M when the core is very large, and somehow not have 
the charge on the core itself directly ejected by the enormous electrostatic forces 
(other than the discharge by the pair production process), the pair production rate 
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would always be more than 26 orders of magnitude smaller than that of a putative 
dyadosphere. In particular, the minimum dyadosphere value would be more than 
1.5 x 10 26 times larger. 

If one wants to have an explicit formula for the pair production rate, this is one 
place where it helps to have the more sophisticated estimate z(l) of Eq. (J9*3|) for 
z(l) rather than just the crude estimate of Eq. (J2*7|) . or the crude estimate z(l) 
of Eq. (jll2j) . since —z is exponentiated in calculating the rate. For example, for 
fj, — 1, using z* instead of z{l) in Eq. (jll7J) would have led to the estimate of 
6.67244 x 10 -27 , which is too large by about 1.7%, and using z(l) would have led to 
8.83141 x 1CT 27 , which is too large by about 33.6%. On the other hand, using z(l) 
instead of z(l) in Eq. (ITTTI) leads to Eq. flU and the estimate p(l) « 6.61149 x 10" 27 
for /i = 1, which differs from the numerical result by less than 0.003%. 

If, for different values of /i = M/M®, we wish to estimate z(l) and the maximal 
ratio, p(l), of the pair production rate to that of a dyadosphere, we can simplify 
the algorithm leading to Eqs. (|9*4^) and (JHHJ) for z{l) and for p(l) by using the fact 
that at the horizon (v = 1) (at least for M < 3 x 1O 6 M ; see below), Eq. (|9~T|) gives 
X 2 (l) greatly dominating over Xi(l) of Eq. (J90|) . so that 

5 15 5 15 
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2zi(l) 16zi(l) 2 2z(l) 16z(l) 2 ' 
Then one gets the simpler explicit approximations for z(l) and for p(l) as 
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57.60469 + 1.932364 ln/i + 0.0011044 In 2 //, 
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p(l) = ^e^exp 
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6.72439 x 10~ 2 V 2 exp 
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16/ 



6.61152 x 10~ 2 V 2 (1 + 0.00055366 lnp - 0.00001764 In 2 p). 



(1: 



Remember from Eq. (JZHJ) that L = In (A/j 2 ) m 65.69511 + 2 In//, with A = 
q 2 m 2 M 2 /ir w 3.39643251 x 10 28 defined in Eq. ©, and with p = M/M , the 
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collapsing core mass in units of the solar mass M . Therefore, a reasonably good 
approximation for p(l) is ir 2 e n /(Ap 2 ) « 6.72439 x 10~ 27 /i~ 2 , which is always within 
about 1.7% of p(l) given above for p, > 1. This is the same approximation as that 
given by Eq. (}2*9"j) of Section 2.2, before the differential equations were derived and 
solved, so it confirms the remarkable accuracy of that approximation. However, as 
we shall see below, both ^e 77 /(Ap 2 ) and p(l) become poor approximations for the 
actual numerical values of p(l) for M > 3.3 x 1O 6 M , where during the collapse the 
maximum electric field remains too weak to give significant discharge, so that the 
self-regulation of the electric field hardly sets in and all of the approximations above 
that rely upon it become poor. 

By comparing the coefficients of the powers of In p in Eqs. ()119|) and ()12U|) (which 
were obtained from the explicit approximations by evaluating and differentiating 
at In /i = 0) with those of Eqs. ()113j) and ()117|) (which were obtained from the 
numerical solution by the slightly different technique of making a fit at p = 1, 
p = a/10, and p = 10), we can get some idea of much error there is in the polynomial 
expansions (in ln/i) of the simpler explicit approximations of (J119)) and ()120j) for 
the normalized inverse electric field strength z = ttE c /E and the normalized pair- 
production rate p = N ' jM c when the core surface enters the black hole event horizon 



For example, the relative error of the quadratic expansion of Eq. (I119j) for z(l), 
compared with the numerical results for z(l), is only about —2 x 10~ 6 at p = 1. 
As p is increased, it becomes positive and rises to a maximum of about 7 x 10 -4 
at p ~ 3 x 10 6 , and then it decreases to cross zero at p « 3 420 991 and thereafter 
rapidly become more negative. It reaches —0.01 (1% relative error) at p ~ 3 660 314, 
reaches —0.1 (10% relative error) at p w 4119 411, and reaches —0.5 (50% relative 
error, too small by a factor of 2) at // ~ 7515 968. 

On the other hand, the pair production rate relative to the minimal dyadosphere 
rate, p(l), is exponentially sensitive to changes in z(t), so the relative error of the 
series expansion for p in Eq. (|120j) . relative to the numerical value, has much more 
rapid changes with /i for p, > 3 x 10 6 . This relative error starts around —2.5 x 10~ 5 at 
p = 1 and at first gradually becomes more negative, going the furthest below zero, 
to about -1.24 x 10~ 3 , at p w 2 740 000. Then the relative error heads back toward 
zero, crossing zero at p ~ 3 134 777 and thereafter rapidly becoming more and more 
positive. It reaches 0.01 (1% relative error) at p ~ 3 292166, reaches 0.1 (10% 
relative error) at p w 3 457659, reaches 0.5 (50% relative error) at p w 3 578 246, 
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and reaches 1 (100% relative error, too large by a factor of 2) at /x ~ 3 627829, 
all before the quadratic expansion in Eq. (J119)) ever leads to a relative error in its 
estimate of z(l) of more than about 0.7%. 

The relative error of the explicit formula p(l) given by Eq. ()99j) at v = 1, 
compared with the numerical calculation of p(l), increases slightly more rapidly 
with fi for n > 3 x 10 6 , reaching 0.01 at p w 3 283 284, 0.1 at /x « 3 456 590, 0.5 at 
/x « 3 577 962 and 1 at /x « 3 627 649. 

The relative errors of the explicit formulae for the barred or hatted quantities 
thus stay small until /x = M/M & gets into the neighborhood of BC ~ 3 695 647. In 
particular, p(l) stays within 1% of the numerical answer until /x reaches 0.8884 BC. 
The error arises essentially because for such large core masses, even though we have 
been assuming Q = M initially at infinite radius (and then ignoring the electro- 
static repulsion and gravitational effects of the electric field), the electric field never 
gets large enough to lead to significant discharge. Therefore, the discharge pro- 
cess never reaches the self-regulation phase in which the solution for z(v) gets near 
the separatrix. Then it is no longer true that Xz(v) ^> X\{v) for v near unity, so 
X(v) = Xi(v) + Xiiv) is no longer an excellent approximation to X(v). For such 
large /x, with little discharge before the collapsing core falls through the event hori- 
zon, it would be better simply to use Xi(v) as the approximation for X(v) for all 
v. 

With a bit of work, one could modify the explicit formula given above to agree 
well with the numerical solution of the ordinary differential equations for larger 
values of /x, but it would be somewhat pointless, since for such large values of 
/x, if the core collapse really did start with Q = M at core radius R = oo, one 
would need to use the Reissner-Nordstrom geometry and also take into account 
the effect of the electrostatic repulsion. On the other hand, for /x < 3 x 10 6 or 
M < 3 x 1O 6 M , there is significant discharge of the core surface by the time it 
reaches the horizon, so by that time the electrostatic repulsion and the deviation 
from the vacuum Schwarzschild geometry become negligible, and the separatrix 
solution becomes an excellent approximation to the electric field and the discharge 
rate. 

For example, if one uses the numerical calculations (assuming a freely infalling 
core surface in the Schwarzschild geometry), one finds that by the time the core 
surface crosses the horizon, £ = Q/M has dropped to £ = 0.1 for /i pa 350 055, 
to £ = 0.5 for n pa 1818 626, to £ = 0.9 for fj, pa 3 319 075, and to £ = 0.99 for 
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/i ~ 3 737 731. If one takes the value of \i that gives £ = 0.9, so that the core has 
discharged only 10% of its charge, one finds that p(l) has a relative error of only 
0.0149, and p(l) has a relative error of only 0.0162, relative to the numerical solution. 
However, the numerical solution itself probably has large errors in this region of 
high £ = Q/M at the horizon, for then the approximations of the Schwarzschild 
geometry and negligible effect of the electrostatic repulsion on the infall will be 
poor. Therefore, the numerical solution itself is probably only valid for /i <C BC ~ 
3.7 x 10 6 , and for this range of fi, the explicit formulae (at least for the barred and 
hatted quantities, though not for the pair production rates obtained from z) are 
excellent approximations to the numerical results. 

Although the numerical solution would not lead to a good estimate of £(1) = 
Q/M when the core surface enters the event horizon (v = 1) unless £(1) <C 1 there, 
it does give an upper limit on the horizon value £(1) as a function of \i = M/M & . 
For example, for fi = BC ~ 3 695 647, the numerical solution gives £(1) ~ 0.98513 
as the upper bound. Lower values of [i lead to lower upper bounds on £, so one can 
say with confidence that to form a black hole with Q > 0.99 M, one needs to start 
with a core that has a mass M > BCMq « 3.7 x 1O 6 M , and probably actually with 
a rather higher mass (even if one could start the core with Q = M when R ^> M). 
I would be surprised if there were any way that an astrophysical black hole could 
form that is within 1% of the charge limit unless it has more than four million solar 
masses. 

Actually, if one were to have a collapsing core sufficiently massive that it could 
collapse into a black hole without discharging significantly, it is very hard to imagine 
how the positive charge (e.g., protons) could avoid being electrostatically ejected. 
Even if nuclear forces were somehow effective in accomplishing that Herculean feat 
for neutron star cores, it would seem even much more unlikely that one could form 
a neutron-star-like core of very many solar masses, so that nuclear forces on the 
protons could conceivably be effective in overcoming the huge electrostatic repulsion 
if there were a significant charge imbalance. 

Therefore, I would actually be surprised if any black holes of astrophysical masses 
ever form within our universe (or our pocket universe with our values of the masses 
and charges of the electron, proton, and neutron) with values of their charges at all 
near their masses. To put it more concretely, I would predict that no astrophysical 
black hole ever has a detectable change in its geometry given by the energy in its 
macroscopic electric field. In other words, Q/M would always be so far below unity 
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that the metric of an astrophysical Reissner-Nordstrom or Kerr-Newman black hole 
would be indistinguishable from a Schwarzschild or Kerr black hole. 

The most plausible way I see to try to maximize £ = Q/M would be for a cluster 
of galaxies, say of mass Mq and with a small excess charge Qc, to collapse to form 
black holes which then Hawking radiate down to M comparable to Qc- 

But if the electrostatic repulsion for protons from the cluster were greater than 
the gravitational attraction, the excess charge in protons would most probably es- 
cape (since protons in objects much larger than neutron stars, such as in galactic 
clusters, surely cannot be all bound by nuclear forces or other forces stronger than 
the gravitational binding like the ones I have been assuming, for the sake of argu- 
ment, in the calculations above of the discharge of a hypothetical collapsing charged 
stellar core). Therefore, one must have qQc < Mcm p . If we are extremely generous 
and say that superclusters of mass Mc ~ 1O 16 M ~ 10 54 are gravitationally bound, 
despite the presently observed gravitational repulsive effects of the dark energy or 
cosmological constant, then we get an upper limit of Qc < 10 36 in Planck units. 
But this is two orders of magnitude below a solar mass, so when Hawking radiation 
reduces M to get near Qc, pair production would radiate the charge away to keep 
Q always much lower than M. 

The only possibility of getting Q ~ M that I can see for a macroscopic black 
hole (i.e., not one that conceivably might be produced in a high energy collider 
from a small number of elementary particles) is to have the initial value of Q larger 
than roughly Q c ~ BCM & ~ 10 44 , but this would require an initial mass of M > 
qQc/iTi P ~ 10 62 ~ 10 24 M Q , which is greater than the mass of all the galaxies in the 
observable universe. Such a great mass does not seem to be gravitationally bound so 
that it could collapse to form a black hole in the presence of the large-scale repulsive 
gravity from the dark energy or cosmological constant. So I would predict that, 
most probably, macroscopic black holes with Q ~ M will never form within our 
pocket universe. 

3 Energy efficiency of the pair production 

We have calculated that even with idealized conditions of a collapsing stellar core 
initially somehow having Q = M and somehow keeping its excess protons from being 
driven off by the extremely large electrostatic forces (though admittedly smaller than 
the nuclear forces within a nucleus), one cannot get the electric field to become large 
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enough to produce pairs at a rate per four-volume within 26 orders of magnitude 
of the minimal dyadosphere rate. However, we did get astrophysically significant 
pair production during this idealized process, and so one might ask what the energy 
efficiency of this process is, what fraction e of the mass-energy M of the stellar core is 
converted into outgoing positrons, to give them total energy eM. Here we shall show 
that unless M > M , the efficiency is very small, e < 1.86 x 10~ 4 ^M/M & < 1. 

In the calculations above of the upper limits of the electric field and of the pair 
production rate (which occur at the surface of the collapsing charged core just as 
it enters the black hole horizon), we assumed that the initial charge Qo of the core 
(when it started what we conservatively assumed was free fall from radial infinity) 
was the maximum possible, Qo = M. Although this indeed would give the maximum 
possible electric field and pair production rate, those quantities at the black hole 
horizon would be nearly independent of the initial charge Qo, so long as Qo was 
large enough to lead to significant pair production and hence to the self-regulation 
of the electric field during the core collapse. 

We found above that for Q Q = M, the explicit approximation for the normalized 
pair production rate, p(l), stays within 1% of the numerically calculated value for 
/i < 0.8884 BC ~ (8/9)BC, indicating that then the self-regulation is sufficient 
to keep the explicit approximation accurate to within about 1% error. If we had 
allowed Qo to be less than M, the analysis would have had the initial value yo of 
the normalized core charge y = qQ/(47rm 2 M 2 ) = hcqQ/(167r 2 eoG 2 m 2 M 2 ) (where 
q and m are the charge and mass of a positron, and M is the mass of the core) 
at core radius R = oo or v = \jlMjR = be not yo = B / ' ji = BM Q /M but 
rather be y = BM Q Q /M 2 . That is, one would just replace B//J, = q/(4irm 2 M) 
with (B/ fi)(Qo/M) = qQo/(47rm 2 M 2 ) in the formulas above. Then the condition 
Bjn > 9/(8C) would become Q /M > 9M/(SBCM Q ) = [9/(8BC)]fi w 3 x 10~ 7 /i- 
For values of the initial charge-to-mass ratio of the collapsing core greater than this, 
there would be sufficient pair production, discharge, and self-regulation that the 
barred approximations above would be fairly good. 

However, the energy emitted into positrons which escape (the electrons going into 
the core or else into the black hole and hence not escaping) does depend significantly 
upon the initial charge Qo (since of course if Qo is lower, there is less charge to be 
emitted), so for this section I shall include that dependence. It is helpful in doing 
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this to define the quantity 

^ = mV)^- = Amn2My = ^ (121) 
M q B 1 y ! 

which is the charge-to-mass ratio at any point in the spacetime outside the collapsing 
core, say labeled by the radial null coordinates (U,V), since the assumed spherical 
symmetry means that there is no dependence on the angular coordinates (0,<f>). 

The initial value of this charge-to-mass ratio £, at past null and timelike infinity 
outside the core, will be denoted by the constant value £ (with maximum value 
1, as it was assumed to be in the previous sections), the value on and outside the 
collapsing core surface when it begins its hypothetical free fall from radial infinity 
in the assumed Schwarzschild metric (ignoring the electrostatic and gravitational 
effects of the electric field, as I have always been doing to get overly conservative 
upper limits on the electric field and pair production rate): 

( ttS({ - o0 ,V) S ^ = !f. (122) 
If one uses ingoing and outgoing radial null coordinates U and V respectively 



that both take the value v = \j2MjR at the surface of the freely collapsing core 
when it has Schwarzschild radius R (the circumference of the core divided by 2tt), 
then the value of £ at the core surface as a function of v may be denoted as £ s , 

«- = «") s «"- b ) = b^oo- < 123 > 

It is also convenient to define £oo(i>) as the value of the charge-to- mass ratio £ at 
radial infinity out along the outward radial null geodesic from the core surface at 



U = v = ^2M/R: 



i 00 = ioo(v)=Z(v,<x>) = ( ^-. (124) 

If a positron is created very near the core surface, it will initially experience the 
electric field of the charge Q s (v) = M£ s (v), which at the core surface will be 

E (v) _ q s (v) _ Qsivy _ uvy (l25) 

Es[v) -RHv)- ~Jm^~ - (125) 

Then as the positron moves out very nearly along the radial null geodesic with 
U — v — const, and V increasing from V — v at the surface to V = oo at radial 
infinity, and as the positron passes the incoming electrons that are created at greater 
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radii and propagate inward to smaller radii, the charge the positron experiences as 
being at smaller radii will decrease to the asymptotic value Qoo(v) = M£oo(t>), and 
the electric field will approach M^ OQ (t; )/r 2 (v , V), which of course becomes smaller 
and smaller as r(U, V) = r(v, V) becomes larger and larger along the outward path 
of the positron. 

In all of this, we make the assumption, to be verified below for M <C BCM Q , 
that only a tiny fraction of the initial total energy of the core is converted into 
outgoing positrons, so that the mass M that appears in the metric is very nearly a 
constant everywhere over the spacetime outside the surface of the collapsing core, 
giving an approximately Schwarzschild metric. 

In this approximately Schwarzschild metric ([30jh as written in terms of radial 
null coordinates U and V, if a positron is produced at some (U, V) and then is 
accelerated to very high gamma factors essentially along the outward null line U = 
const., then the kinetic energy gained during this acceleration by the electric field 



where the integral is taken along the outward radial null line U = const, along 
which the positron approximately travels from its creation point at (U, V) to radial 
infinity at (U,oo). Almost all of the positrons are created sufficiently deep in the 
electric field that their final kinetic energies, at radial infinity, far exceed their rest 
mass energies, so I shall ignore the latter. Then one may multiply the positron 
production rate by the energy gained by each positron and integrate over all of the 
exterior region where the production is occurring to get the total energy emitted 
during the core collapse as 



(The negative sign comes from the fact that Q is decreasing during the discharge, 
so Q >uv < 0.) 

In principle one could insert the approximation z(v) from Eq. (JH%|) into Eq. 
(JBUj) to get an explicit approximation for the distribution of charge Q(U, V) in the 
region outside the collapsing core where the dominant pair-production occurs, and 
then insert these results into Eqs. (|126j) and (J127)) to do numerical integrations over 
(U, V) to find the total energy eM emitted in positrons. For even more precision, 



Q(U,V)/r 2 (U,V) is 




(126) 




(127) 
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one could directly solve the pair-production partial differential Eq. (|37|) numerically 
and to determine the integrands of Eqs. 1)126)1 and (J127j) . but this shall be left as 
an exercise for the reader. However, here I shall resort to approximations, which, at 
the cost of some precision, have the advantage of giving explicit results. 

Just as the derivation of the approximate ordinary differential Eq. (|77j) for the 
discharge of the core from the partial differential Eq. (J37|) relied upon the fact that 
z" 1 <C 1, so I shall use this same fact here to make an approximation valid only to 
zeroth and first order in z~ x . This fact also implies that, for v not too close to unity, 
the almost all of the pair production occurs close to the collapsing core surface, near 
r = R, and this implies that £(U,V) = Q(U,V)/M is nearly constant along the 
outgoing radial null lines, equalling both £ s (i>) and £oo(v). 

An exception to this is for 1 — v <C 1, when the core surface is just about to 
enter the black hole horizon with finite nonzero charge-to-mass ratio £(1,1), but 
when the pair production outside the black hole reduces the charge-to-mass ratio 
at radial infinity, £oo( v ) — £(^5°°); to zero in the limit that U = v — > 1, which is 
the limit of U = t — r* — > oo at future null infinity, assuming for this argument 
that the black hole lasts forever and does not evaporate away by Hawking radiation. 
That is, for a positively charged black hole that lasts forever and does not have any 
incoming charge other than the electrons produced by its electric field, eventually 
the charge will be seen to decay to approach arbitrarily close to zero, because of the 
continuing pair production outside the black hole and resulting flow of positrons to 
radial infinity. 

(Strictly speaking, even this is not quite right, since the outgoing positrons will 
not have quite the speed of light and hence will actually go to future timelike in- 
finity rather than to future null infinity, so that along future null infinity itself, the 
charge will always remain at its initial value Qq that it has along past null infin- 
ity and spacelike infinity. However, at any large but finite function of the 
Schwarzschild time coordinate t, which is proper time for a static observer at radial 
infinity, the charge will asymptotically decrease to zero as the external electric field 
of the initially- charged black hole produces charged pairs, with the electrons falling 
down the hole to neutralize it and the positrons carrying off the charge to radial 
infinity, even though actually at future timelike infinity rather than strictly being 
at future null infinity. Since under the conditions that the charge self-regulation is 
significant, the positrons are accelerated to enormous gamma factors by the elec- 
tric field, for this paper I am employing the approximation that the positrons go 
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out at the speed of light and so am ignoring the distinction between their outgoing 
worldlines and outgoing null geodesies.) 

For the formulas below, I shall assume that £ > [9/(85(7)]// ~ 3 x lCT'/i, so 
that the self-regulation of the electric field is significant even before the collapsing 
charged core gets near the black hole horizon, so that most of the decrease of £00 (^) 
from £0 at v = to at v = 1 occurs for v not too close to 0. Then the self-regulation 
of the charge by the pair production process gives ^ (v) <C £0 before v gets so close 
to unity that the charge-to-mass ratio at radial infinity along the outgoing nearly- 
null positron worldlines, £°o(f), drops significantly below the charge-to-mass ratio 
£ s (f) at the core surface at the start of these outgoing positron worldlines. 

In the limit that this approximation is exact, that is, to zeroth order in z -1 , one 
gets that S(U, V) = qQ s /R = \q^ s {v)v 2 , This would give 
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However, the charge Q and the charge-to-mass ratio £ both actually decrease a 
bit as one moves outward along each outward null geodesic U = v = const.. If one 
uses Eq. (I60J) (with P replaced by P that is the improved value for v <C 1, though 
I shall remind the reader that I do not know the improved value for v ~ 1 and so 
there am likely to have errors of the order of v/z) for the charge distribution outside 
the collapsing charged core surface, one gets that 
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+ P(v) + 1) 2 

where here X = (V - U)/2 = (V - v)/2 « M(r - R)/R 2 is the spatial radial 
coordinate X given by Eq. (|47jl. not the X = (1 — Y)/Y of Eqs. (jl()2j) and onward, 
and where Eq. (jlOO)) implies that the P defined in Eq. (176)1 may also be written in 



terms of v and of Y(v) 
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defined in Eq. (|100|) . as 



One can see from this that 
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is thus a bit less than £o(^)> by an amount that is 0{z^ x ) multiplied by a logarithm. 
When £ > [9/(85(7)]// ~3x 10~ 7 yU, then for sufficiently large v, (, s (v) will decrease 
significantly below £o, and in this large v region, the electric field self-regulation will 
lead to Y pa 1 and hence to 

U(V) « -jr^ f 1 7~T l n T^—) ■ (132) 

Bv^zyvj \ z(v) 1 — v J 

This equation would imply that £oo(^) would become negative for t> > 1 — e~ z( - v ^ 2 , 
which is for v very near unity (very near the point where the collapsing charged core 
surface enters the black hole horizon at v — 1), but actually the approximation is 
breaking down for v this close to unity. As discussed above, one would instead expect 
that £oo{v) should asymptotically approach as v — > 1, as then there is an infinite 
amount of time in the region external to the black hole for the pair production 
to carry away all of the charge (under the assumption that the black hole persists 
forever and does not evaporate by Hawking radiation, other than the charge emission 
that is of course included in the full Hawking emission formula). Therefore, we shall 
take Eqs. (|129j) - (jl32|) to be approximately valid only for — In (1 — v) <C z(v)/2. 

If we define 

P(v) , \ 

a(v) = y ^ w v, (133) 



where the ~ applies only for the large-f region where £ s (f) drops significantly and 
Y pa 1 (the region where the pair production is significant), then we can calculate 
that for fixed U = v , the mean value of the spacelike radial coordinate X(U, V) = 
X(v, U) at which positrons are created along the outward radial null geodesic is 
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W ~ 4(1 - v)z{v) In [1/(1 - a)} ^ n 2 ' (134) 
Then r pa R(l + RX/M) = R(l + 2X/v 2 ) gives 

(1\ ~ i + F(v) = l + V — 

K R' {) 2(1 - v)z(v) In [1/(1 - a)} ^ n 2 

1 00 })" 

1 + 2(1 - v)z(v) In [1/(1 - v)} § ^ ' (135) 

where again the last ~ applies only for the large-f region where £ s (v ) drops signif- 
icantly and Y ~ 1, so that a(v) ~ v. For d C 1 and yet still within the region 
where £ s (i>) has dropped significantly below £ > one gets (r — R) pa R/(2z), just as 
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was found in Section 2.3 for getting the approximate estimate of the discharge rate 
before solving any differential equations. The main correction for somewhat larger v 
is to divide this result by 1 — v, making the mean radius of the positron production 
somewhat larger. This approximation would imply that the mean radius would go to 
infinity as v — > 1, but this is where the approximation is breaking down. However, 
for £ > [9/(8BC)]fi ~ 3 x lCT 7 ^ so that significant self-regulation occurs, most 
of the pair production occurs for v sufficiently below unity that the approximation 
does remain valid, so for calculating the efficiency e in such cases, one does not need 
to worry about the breakdown of the approximation for 1 — v -C 1. 

Now when 1—v ^> 1/ (2z) ~ 0.01, an estimate for the average energy (for given v) 
gained by a positron in being electrostatically accelerated out to radial infinity that 
is better than the zeroth-order (in z~ l ) estimate qQ s /R = qM^ s / R is the first-order 
estimate 

mn>w i^„2^g, (136) 

because most of the work done by the electric field in accelerating the positron 
from its average creation position at r = (r) »* R{1 + F) to r = oo occurs for 
r - i? > (r-R) where Q(U, V) « Q^v) = M^(v). 

One can now insert this expression into Eq. (|127|) to get an approximate one- 
dimensional integral for the energy efficiency of the pair-production process (depend- 
ing on the core mass M .!/•// and on its initial charge Q = M£ = M © /i£o) : 

q Jo dv Jo (r) 

1 f 1 v*Uv) dUv)_ dv m 



2 Jo 1 + F(v) dv 

It was this integral that was evaluated numerically from the numerical solution of 
Eq. (|77jl. using the approximate expressions for ^oo(f) given by Eq. ()131j) and for 
l + F(v) given by Eq. $35$ . 

However, before Eq. (|137|) was integrated numerically, an explicit approximate 
solution was found. For this it appeared too complicated to use the approximation 
z of Eq. ([SI)) , so I used instead an alternate version of z{y) that was given by Eq. 
(j!12jl . This alternate simple approximation is that z(v) is given approximately by 
z(v) that is the solution to the equation 

zV = V^" 5 + B^^v-'e 3 - 1 ^ 1 -- 4 . (138) 
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Here, unlike in Eq. ()112|) . I have allowed the possibility that the initial charge-to- 
mass ratio £o is less than unity by replacing Bj\x there by -B£o/a* = qQ / (Airni 2 M 2 ) . 

In this expression, the second term on the right dominates when v is sufficiently 
small that the pair production is small. Indeed, if one dropped the first term as 
negligible compared to the second, this would give z ~ B^fi^v ~ 4 , so then the 
normalized charge would be y = v~ 4 z~ 1 « v~ 4 z~ l m B£q/[i = yo = const.. There- 
fore, this regime would give £ s (v) ~ const, and hence a negligible contribution to 
the integral (j!37|) for the efficiency e. That is, the integral for the efficiency is domi- 
nated by the values of v for which the first term on the right hand side of Eq. ()138)) 
dominates, the region of significant charge self-regulation. 

If we define 

5 5 3 5 

L = L + -ln£o = ln,4 + -lnB + -ln/i + -ln£o 

= -- In (2 W 7c 9 m 2 M 2 q- 13 Q 5 ) « 79.01593 + 0.75 In // + 1.25 In £ Q , (139) 

then the two terms on the right hand side of Eq. (J138)) are equal at the transition 
value of v that I shall call vt = Vt([i, £o) and which is given by the following equation 
in terms of the related quantity z t = z t (n,£, ) that is defined by the next equation 
after that: 

/ \ 1/4 

' (140) 

3 

zt + -In z t = L . (141) 

Under the condition that significant self-regulation of the charge occurs, that is, 
for £ > [9/ (8BC)]fi ~ 3 x 10 -7 /i, a sufficiently good approximation to the solution 
of Eq. is 

3 - Q „ „ 27- 
zt « L - - In A) + -L, 1 lnL + — L ~ 2 lnL (lnL - 2), (142) 

which one can then insert into Eq. ()140|) to get v t . Two rougher but more self- 
contained expressions for z t are 

zt~ - hn(2 10 7i 9 m 2 M 2 q- 13 Q 5 ) 



-iln(2 10 7r 9 m 2 M 2 g- 13 Q 5 ) 



75.77 + 0.7426 In /i + 1.238 lnf . (143) 
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One can see that when v = Vt, both terms of the right hand side of Eq. ()138|) 
are equal to z 2 e Zt , so one gets 

z(v t ) + 2 In z(v t ) = z t + 2 In z t + In 2. (144) 

Therefore, the transition value of z defined in this way is just a bit larger than z t , so 
Zt is a fairly good estimate of the value of z where the self-regulation of charge starts 
to become important as the surface of the charged core freely collapses inward from 
the idealized initial conditions of being at rest at radial infinity and then follows 
radial geodesies inward in the assumed exterior Schwarzschild metric of mass M. 

For ii — 1 (M = M Q , one solar mass) and £ — 1 (Qo — M, initial charge equal 
to the mass so that then one would have the extreme Reissner-Nordstrom exterior 
metric, though the infall is always taken to occur in the Schwarzschild metric with 
no effect from the electrostatic forces), one gets from Eqs. ()142j) and ()140|) that 
z t (fi,^ ), v t ([i,!;o), and i?i(/i, £ ) take the values 



z m = z t (l, 1) « 75.77015, (145) 
v t = v t (l, 1) « 0.02360985, (146) 

% = = -i- « i793 - 964 - ( i47 ) 

2M M v 2 n y > 

In fact, for general \i = M/M Q and £ = Qo/M, the critical radius where the 
self-regulation starts becoming significant is 



Rt = 2M/vf w 5298.024^0^/75.77015 km ~ 5298 ^- 5049 ^- 5082 km. (148) 

Thus if one starts with a hypothetical charged core collapsing freely from a very 
large radius into a black hole with initial charge Qq ~ M, the self-regulation of the 
charge will start to become important when the the core gets to a radius of the 
order of the radius of the earth, and that is when there will start to be significant 
pairs being produced. The proper time left during the free-fall collapse after the 
core surface crosses this radius is then 

At = -M/vf = -M Q (B 3 ^z*) 1/4 ~ 0.499 fi - 257 ^ 762 seconds. (149) 
3 3 

Therefore, if Qo ~ M, the burst of positrons that will be emitted as the core collapses 
from R = R t to R = 2M lasts a time that is of the order of one second. 
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For v > Vt, where the pair production is significant, the first term on the right 
hand side of Eq. ()138|) dominates, and one gets z 2 e z ~ Afi 2 /v 5 . To evaluate the 
integral in Eq. ()137|) . it is simplest to express z as a function of the surface charge- 
to- mass ratio £ s , which leads to the equation 

z+- A hxzPaL = L Q + ~hx^- : (150) 
4 4 4o 

with the approximate solution 

3 ~ 9 - - 27 - - 

z ss L InL H LT l InL + L~ 2 InLflnL — 2) 

4 16 128 v ' 

~ Zt fl + -A_i n |lV (151) 



4zt + 3 6 



o, 



Then when one evaluates the integral in Eq. ()137|) to first order in z t _1 , one gets 
that the efficiency of the conversion of the core mass M into outgoing positrons of 
energy eM is, from the explicit approximation, 



12^ 5z t 1 4 j 

1 «*2 2 1 6 ^ 



pa ^^o)-3«(l ^ ^ 
« 0.000185534 ^/^o 74 (1 - 0.004901 ln/i - 0.008168 ln£ ) 
« 0.0001855 /i°- 4951 e - 7418 . (152) 

By comparison, a numerical evaluation of the integral in Eq. (J137J) from the 
numerical solution of Eq. (|77jl for £o = 1 an d for \i — 1 and for /i = 10 and then fit 
linearly to /i _1 / 2 e(/i, 1) versus ln/i gave 

e(/i, 1) w 0.000185467 /x 1/2 (l - 0.004980 ln/i). (153) 

This shows that over a range of /i of the order of unity, the explicit approximate 
expression seems to agree to about 4 decimal places with the numerical calculation, 
so it is apparently quite accurate. However, I have noted above that even my 
numerical calculation is based upon certain approximations to reduce the partial 
differential equation (|3*Tj) to the ordinary differential equation (|77jl . and I would 
expect that the relative error is of the order of v/z. For /i = 1 and £o = 1, most 
of the positrons are produced around v ~ v t n ~ 0.024 and z ~ z t \\ ~ 76, giving 
Vtn/ztll pa 3 x 10~ 4 as a crude estimate of the relative error. Therefore, I might 
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guess that the results for the energy efficiency above are accurate to about 3-4 
decimal places and can be written to roughly this precision as 

'M\ "' 5 /-«o- ' 742 



e « 0.0001855 ^-j . (154) 

To show the main dependence of the efficiency upon the mass m and charge — q 
of the electron and upon the mass M and initial charge Qo of the charged core that 
collapses freely from rest at infinity in the assumed external Schwarzschild metric, 
let us take the first term on the right hand side of Eq. (jl42j) or Eq. (j!43|) as the 
crude approximation Lq for z t and then insert this into Eq. ()152j) . with the 0{z^ 1 ) 
corrections omitted, to get, with conventional units restored, 



m 

e ~ 



9° / 15g x 

3M^ceogln[(2 2 % 18 ^ 7 c 7 e9G 2 m 2 g~ 13 M 2 Qo 5 )- 1 /4]- v 1 

For M = M Q « 1.98844 x 10 30 kg and Q = y/Aire GM e « 1.71353 x 10 20 coulomb, 
this crude estimate gives e ~ 0.000182, within 2% of the efficiency 0.0001855 ob- 
tained numerically for a solar- mass core with maximum initial charge, £o — Qo/M = 
1 in Planck units. 

Another way to express the efficiency is to note that Eq. ()152|) may be written 

as 



m 



4vrm 2 M / 1 Qv t \ 1/2 3/4 



3m p J \ yfazt V 12z t hz t 



(0.0001815) (1.0235) (0.9985) Ai - 4951 ^- 7418 



(0.0001815) (1.0220) /i°- 4951 e ( - 7418 



o 

7418 
i 



w 0.0001855 /i u - 4y&i e o ams , (156) 

where (l/3)m/m p gives the factor of 0.0001815, the square root (when z t is evaluated 
to give 75.770 for Q = M = M ) gives the factor of 1.0235 that is accidentally 
close to unity, and the truncated series with the 0(z t -1 ) corrections gives the factor 
of 0.9985. Therefore, the main cause for making the very conservative upper limit 
on the efficiency e much smaller than unity is the small ratio of the mass of the 
electron to the mass of the proton, and it turned out coincidentally that for a solar 
mass core, the maximum idealized efficiency of converting maximally charged core 
energy into outgoing positron energy is within about 2.2% of one-third this small 
mass ratio. 
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If one multiplies the efficiency e from Eq. ()154|) by the mass-energy M = Mc 2 
of the core, one gets that the total energy emitted in positrons is 

f M\ 1,495 fd \ °- 742 
eM « 3.315 x 10 50 f — J [^j ergs, (157) 

which for a core of mass M = m~ 2 ~ 1.85327M Q and of maximal initial charge 
(Qo = M = M) is very roughly mj (3?7ip) = hem/ (3Grrip) ~ 6 x 10 50 ergs. For stellar 
mass cores, this is 2-3 orders of magnitude smaller than the energy of gamma-ray 
bursts |H ES] , essentially because the efficiency is 3-4 orders of magnitude less than 
unity. 

Again, I should emphasize that all of these estimates and calculations give only 
very conservative upper limits on the efficiency and energy emitted, since they all 
assume that somehow one can hold the charge onto the surface of the collapsing core 
even when it is as large as the transition radius Rt, which for £o — Qo/M > 10 -4 
would be significantly greater than the size of a neutron star. For such a core the 
electrostatic forces of repulsion on the excess protons on its surface would be more 
than 14 orders of magnitude greater than the gravitational attraction of the core for 
these protons, and it is hard to imagine that any forces sufficiently powerful to hold 
in the protons (such as nuclear forces at nuclear distances) could be effective in any 
star larger than a neutron star. And if one does take £o ~ 10 -4 so that R t ~ 50 km, 
somewhat larger than a neutron star size, then the upper limit on the efficiency of 
the pair production process would be only of the order of 2 x 10 -7 , which is far too 
small to give a viable model for gamma ray bursts. 



4 Annihilation probabilities 

We have taken an idealized model in which originally there is no matter outside the 
positively charged collapsing stellar core, but only an electric field that produces 
electron-positron pairs there. We have assumed that the electrons produced propa- 
gate freely to the core surface, and that the positrons produced propagate freely to 
radial infinity (both being accelerated by the electric field, of course). Now we shall 
confirm that indeed the annihilation probabilities are very small for these electrons 
and positrons, so that indeed they propagate essentially freely. 

The cross section for an electron to annihilate with a positron that has a large 
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gamma-factor 7 ^> 1 in the frame of the electron is 

<7«-^fV(2 7 )-l]. (158) 
777/7 

From this, let us calculate the probability P that a positron created very near the 
collapsing stellar core surface annihilates with any one of the incoming electrons as 



the positron propagates out to radial infinity, as a function of v — J2M/R. Since 
the annihilation probability will turn out to be extremely small, it is sufficient to 
confirm that by taking the Newtonian limit (strictly speaking valid only for t)<l, 
but here applied to all v < 1) for the core collapse and exterior geometry. 

Since the electrons and positrons are very rapidly accelerated to near the speed 
of light, in the lab frame (the center of mass frame of the collapsing core) they will 
have a relative velocity of very nearly 2c = 2, so the annihilation probability (if 
small) would be 

POO 

P^2 Fadr (159) 



R 

if F were the inward number flux of mono-energetic electrons. 

However, the cross section depends on the energies of the outgoing positron 
and the incoming electrons, so one needs to divide up the flux of electrons into 
partial fluxes according to their energies. In particular, consider electrons that are 
produced between r e and r e + dr e and then propagate inward to give the partial flux 
(dF/dr e )dr e , where for r e — R <C R 

dF_ m 4 e ~ w ^ _ g 2 Q{t) 2 (_ vrmV|\ 

m 4 e~ z f 2z 



cxp 



4tt z 2 

Then the annihilation probability may be written as 



(160) 



poo ^00 

P^2 dr dr e —a. (161) 
Jr Jr dr e 

At a possible annihilation point r between the core radius R (where the positron 
is assumed to be produced, to maximize the annihilation probability P) and the 
electron production radius r e , the positron and electrons will have accelerated to 
the following gamma factors, respectively, in the lab frame: 

qQ ( 1 1\ qQ(r — R) Tim, /„ „„s 

7+^ — (-=-- ~ „2 J *—(r-R), 162 
m \R r) mR A z 
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qQ (1 l\ qQ{r e -r) irni 

7- « — ~ w 2 — « r e -r , 163 

m \r r e J mR A z 

where the final two expressions on the right hand side of each of these approximate 

equations are valid for r e — R <C R, which is where almost all of the electrons are 

produced. For 7 + ^> 1 and 7_ 3> 1, the relative gamma-factor between the electron 

and the positron, which goes into the cross section formula ()158|) . is 

7 » 27+7_. (164) 

Now we may see that the exponent in the final exponential of Eq. (|16U|) is 
—5(7+ + 7-), where 

$ = (165) 
Then in Eq. f!161|) we may also use r w R + (z/7rm)7 + , r e « r + (z/Trm)'-/-, 
qQ = 4irm 2 R 2 /z, and i? = 2M/v 2 , and also insert 7 from Eq. ()164|) into the cross 
section formula (|158|) to obtain 

P ~ t£L [°°tl± /° ^(ln 7+ + ln7- + 21n2-l)e- 5 ^-) 



4 —2 

g e 


/■°°d7+ /• 


47T 2 J 


1 7+ A 


4 —2 

g e 


, ixmM 

2 _ 2 — 


4vr 2 



+ (2 In 2 - l)ln 



2 nrnM 



V 2 Z 2 



(166) 



When the self-regulation of the charge is effective, the first term on the right 
hand side of Eq. ()138j) dominates, so one has 

*v « 4£ s «!^. (167) 



Then 



_ q 2 v 5 z 2 



nmM „ ^ N , 9 irmM 



In 3 — — + (2 In 2-1) In 2 



47rm 2 M 2 L f 2 -2 2 ' v ' f 2 z 2 J ^ 

The maximum probability this formula gives is when M is minimized and v is 

maximized, so let us set M = M Q as a very conservative minimum value of the mass 

of a collapsing core that can form a black hole, and also take v = 1, its maximum 

value outside the black hole. Then the numerical calculations described above gave 

z ~ 57.6048, which leads to an annihilation probability of only 

P « 3.23 x 1(T 27 , (169) 

which is utterly negligible. Therefore, one may completely neglect the annihilation 
of the electrons and positrons produced by the electric field of a collapsing stellar 
core. 
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The main reason that the annihilation probability is so small is that when the 
positron and the electron have each accelerated a radial distance in the electric 
field that is much greater than the distance m/(qE) = z/(irm) needed to give each 
particle a kinetic energy equal to its rest mass, so that the positron and electron 
energies are E + ^> m and £L ^> m, the annihilation cross section given by (j!58j) . 

7TQ 4 / AE + E^ \ 

1 (170) 



2E + E_ V m 2 

drops roughly inversely with each acceleration distance. Therefore, the integration 
over these acceleration distances gives a result only logarithmic (with the leading 
term going as the cube of the logarithm) in the effective total available acceleration 
distance, which for the region outside a charged sphere of radius R is roughly Rj (2z). 
Even for small z, so that the e~ z factor in Eq. (j!66|) does not give significant 
damping, one needs a very large available acceleration distance (or very small 5) 
in order for the cube of its logarithm to compensate for the q A = a 2 factor in Eq. 

(HEED- 



5 Conclusions 

It does not seem to be at all nearly possible to have astrophysical dyadospheres 
(electric fields larger than the critical value for Schwinger pair production, over 
macroscopic regions much larger than the regions of high fields that might conceiv- 
ably be produced by individual heavy nuclear collisions). If, as is most plausible, 
charge carriers like protons are bound to an astrophysical object, such as a star 
or stellar core, primarily by gravitational forces, then the electric field cannot get 
within 13 orders of magnitude of the minimal dyadosphere values. (First the excess 
charges will simply be ejected by the electrostatic repulsion when that exceeds the 
gravitational attraction. Then pair production rates from the macroscopic electric 
field, as opposed to that from collisions of individual particles, will be trillions of 
orders of magnitude below dyadosphere values and so will be completely negligible.) 

Even in the implausible scenario in which the excess charge carriers are bound by 
nuclear forces to a collapsing stellar core, I have shown here in a simple spherically 
symmetric model that the electric field has a very conservative maximum value that 
is more than a factor of 18 below the minimal dyadosphere value. Because the pair 
production rate is essentially exponential in the negative inverse of the electric field, 
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the upper limit of the pair production rate in even this implausible scenario is more 
than 26 orders of magnitude below the minimal dyadosphere values. 

The idealized implausible scenario considered here, to give this very conservative 
upper bound on the pair production rate, had the maximal amount of charge some- 
how bound to the surface of an idealized stellar core with maximal initial charge 
that undergoes free fall collapse from radial infinity in the Schwarzschild metric 
(conservatively ignoring the fact that an actual astrophysical collapse would start at 
finite radius and fall in slower, giving more time for discharge, and the fact that if 
the initial charge were maximal, the strong electric field would modify the geometry 
and also give electrostatic repulsion of the core surface, both of which would also 
slow down the collapse and lead to greater discharge and smaller electric fields). 
This scenario led to the maximal electric field (occurring when the freely collapsing 
core enters the event horizon of the black hole that would form) being less than 
5.5% of dyadosphere values. Although the pair production rate is always less than 
10~ 26 that of dyadosphere values, in this implausible scenario of having no other 
mechanism of discharging the core, there would be enough pair production to keep 
the electric fields always more than 18 times smaller than dyadosphere values. 

Although astrophysical dyadospheres do not form, in this idealized implausible 
scenario there is significant pair production (though at astrophysical rates that are 
much, much lower than dyadosphere rates), and it may be of interest to calculate 
what fraction of the total mass-energy M of the collapsing stellar core would be 
converted into pairs. Here it was found that this efficiency, even under the highly 
idealized conditions of having maximal initial charge at such large radii that it seems 
inconceivable that the charge carriers could be sufficiently bound to such objects so 
much larger than neutron stars, was always much less than unity for collapsing 
objects with much less mass than three million solar masses: the efficiency was very 
conservatively bounded by 2 x 10~ 4 ^/ M/M & . Therefore, even these idealized charged 
collapsing objects, unless they were enormously more massive than the sun, would 
not produce enough energy in outgoing charged particles to be consistent with the 
observed gamma ray bursts. 

It would of course be of interest to calculate the upper limits on the electric 
field and on the maximum pair production rate for models in which one relaxed the 
spherical symmetry. Although I would readily admit that I do not have a rigorous 
proof that the pair production rate cannot be higher then, the general arguments of 
Subsection 2.3 strongly suggest that it would be very surprising if it could be larger 
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by more than about one order of magnitude. Therefore, the spherically symmetric 
model analyzed here compel me to conjecture that the maximal pair production 
rates achievable by macroscopic astrophysical electric fields are more than 25 orders 
of magnitude below that of hypothetical dyadosphere values. 

In conclusion, dyadospheres almost certainly cannot form astrophysically, and 
the much weaker pair production rates that might occur, under highly idealized and 
implausible scenarios, do not seem sufficient for giving viable models of gamma ray 
bursts. For a shorter version of this work for a conference proceedings, see [35J. 
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